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A  computer  program  was  written  for  use  in  sizing  gas  turbine  counter- 
flow  regenerators.   The  program  simulates  the  semi-open  gas  turbine  re- 
generative cycle  under  constant  load  conditions.   Inputs  to  the  program 
include  compression  ratio,  turbine  inlet  temperature,  component  efficien- 
cies, horsepower  output  and  regenerator  core  geometry.   External  re- 
generator core  dimensions  are  varied  and  various  outputs  are  computed 
such  as  thermal  efficiency,  mass  flow  rate,  total  relative  pressure 
drop,  regenerator  effectiveness  and  temperatures  and  pressures  at  all 
key  points  in  the  cycle. 
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1.    Introduction, 

The  design  of  a  heat  exchanger  involves  a  consideration  of  both 
the  heat  transfer  rates  between  the  fluids  and  the  mechanical  pumping 
power  expended  to  overcome  fluid  friction  and  move  the  fluids  through 
the  heat  exchanger.   Kays  and  London  state  that  for  a  heat  exchanger 
operating  with  high  density  fluids,  the  friction  power  expenditure  is 
generally  small  relative  to  the  heat  transfer  rate,  with  the  result 
that  friction  power  expenditure  is  seldom  of  controlling  influence  (1). 
However,  for  low-density  fluids,  such  as  the  gases  flowing  through  a 
gas  turbine  regenerator,  it  is  very  easy  to  expend  as  much  mechanical 
energy  in  overcoming  friction  as  is  transferred  as  heat.   The  regenera- 
tor designer  must  therefore  establish  the  proper  balance  between  heat 
transfer  rate  and  pressure  drop. 

Two  advantages  of  the  gas  turbine  are  compactness  and  low  weight 
whereas  heat  exchangers  are  inherently  rather  bulky  and  heavy.  Weight 
and  size  considerations  must  therefore  be  taken  into  account  along  with 
heat  transfer  rate  and  pressure  drop  balance  if  regenerators  are  to 
serve  a  useful  purpose.   Since  the  prime  reason  for  regenerators  is 
the  improvement  of  overall  thermal  efficiency,  it  appeared  that  a  fruit- 
ful approach  to  regenerator  sizing  and  core  selection  would  be  considera- 
tion of  the  gas  turbine  cycle  as  a  whole  rather  than  the  regenerator 
separately.  Adjustments  in  internal  regenerator  core  geometry  and  ex- 
ternal dimensions  are  then  reflected  directly  in  terms  of  changes  in 
thermal  efficiency  rather  than  indirectly  in  terms  of  changes  in  pres- 
sure drop  and  heat  transfer  rate.   A  computer  program  was  written  for 
this  purpose  which  essentially  simulates  the  semi-open  gas  turbine  re- 
generative cycle  under  constant  load  conditions.   Inputs  to  the  program 
are  compression  ratio,  turbine  inlet  temperature,  component  efficiencies, 
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horsepower  output,  regenerator  core  geometry,  etc.   External  regenerator 
core  dimensions  are  then  varied  and  various  outputs  are  computed  such  as 
thermal  efficiency,  mass  flow  rate,  total  relative  pressure  drop,  re- 
generator effectiveness  and  temperatures  and  pressures  at  all  key  points 
in  the  cycle. 

Since  the  program  is  based  on  fundamental  relationships  from  thermo- 
dynamics, heat  transfer  and  fluid  mechanics,  sections  2,  4  and  5  are  a 
review  of  the  basic  theory  utilized. 
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2.  Thermodynamics  of  Semi-Open  Regenerative  Cycle 
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2.1  Description  of  the  Actual  and  Ideal  Cycles. 

Both  the  ideal  and  actual  regenerative  cycles  are  shown  in  Figure 
1.   The  state  points  indicating  the  ideal  cycle  are  subscripted  "s"; 
the  nonsubscripted  state  points  indicate  the  actual  cycle. 

First  describing  the  ideal  cycle,  air  continuously  enters  the  com- 
pressor in  condition  "Is".   It  is  compressed  isentropically  to  state  "2s" 
then  enters  the  cold  side  of  the  regenerator.  Theoretically,  if  the  re- 
generator is  large  enough  and  the  flow  is  slow  enough,  the  air  from  the 
compressor  can  be  heated  reversibly  to  temperature  "5s"  at  state  "3s". 
It  next  enters  the  combustor,  some  of  it  going  around  the  outside  of 
the  combustor  chamber  proper  and  some  furnishing  oxygen  for  burning  fuel 
which  is  continuously  injected.  Because  of  the  temperature  rise,  the 
gases  expand  and  enter  the  turbine  in  state  "4s".  After  isentropic  ex- 
pansion through  the  turbine,  the  hot  exhaust  gases  enter  the  hot  side  of 
the  regenerator.  Again,  if  the  regenerator  is  large  enough  and  the  flow 
slow  enough  the  gases  are  cooled  reversibly  to  temperature  "2s"  at  state 
"6s"  and  exhausted. 

The  ideal  cycle  can  be  described  as  follows: 

&<-£,)  r„-s„  (2-2) 

In  the  actual  cycle,  air  continuously  enters  the  compressor  as  in 
the  ideal  cycle;  condition  "1"  is  identical  to  condition  "Is". 

R-R,  (2.4) 
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T>T7s 

(2.5) 

It  is  then  compressed  in  an  irreversible  adiabatie  process  to  state  "2" 
arriving  at  the  same  pressure  as  in  the  ideal  case  but  with  increased 
entropy  due  to  fluid  friction. 


Vi  *  Rs-  (2.6) 


There  will  of  course  be  a  pressure  drop  as  the  air  flows  through  the 
regenerator  and  the  air  is  heated  only  to  the  temperature  at  state  "3" 
not  "3s"  because  less  heat  is  actually  transferred  than  could  be  trans- 
ferred reversibly.  The  air  and  now  gases  undergo  another  pressure  drop 
in  flowing  through  the  combustorj  however  they  arrive  at  the  turbine 
inlet  at  the  same  temperature  as  in  the  ideal  case. 

To  continue  the  analysis  it  is  necessary  to  look  ahead  to  state  point 
"6".  The  pressure  at  this  point  can  be  considered  constant  because  it 
is  assumed  that  the  turbine  is  exhausting  into  a  large  volume,  most 
likely  the  atmosphere. 


R-R: 


(2.8) 


There  will  of  course  be  a  pressure  drop  through  the  hot  side  of  the  re- 
generator which  indicates  that  P_  must  be  higher  than  P_   if  P,  is  to 
remain  the  same  as  in  the  ideal  case.   The  gases  will  also  experience 
another  increase  in  entropy  as  they  pass  through  the  turbine.   Finally, 
since  less  heat  is  actually  transferred  from  the  hot  side  of  the  re- 
generator than  could  be  transferred  reversibly,  T,  is  higher  than  T,  . 
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2.2  Compressor. 

If  the  compressor  is  considered  a  steady  flow  machine  with  kinetic 
energy  change  equal  to  zero,  the  works  of  the  fluid  for  adiabatic  flow 


are: 


ideal:       WCs  »  Klf-  his 
actual:      Wc    -  h.\-fr-L  *  /us  - Im. 


(2.9) 
(2.10) 


The  foregoing  equations  may  be  obtained  from  the  steady  flow  energy 
equation  (2).   The  compressor  efficiency  is  then: 

Jl    _   IDEAL  WA><4»< 


—  ^'5  ~  ^-IS 
ACTUAL  FLUID  WOUK       '     ^(J  __  ^ 


2.3  Combustion  Process. 


k3 


(2.11) 


Figure  2  Energy  Diagram  for  Combustor 

The  energy  diagram  shown  above  for  a  steady  flow  combustor,  is  on 

the  basis  of  one  lb  air  entering  and  k^-A  lb  fuel/lb  air.  The 

m  /mm 

entering  molecular  internal  energy  and   flow  work  of  the  air  is  h_  and 
of  the  combustible  is    Yklo<  A-£     •     The  stored  chemical  energy  in  the  fuel 
q/>i8  measured  from  the  zero  datum  of  enthalpy.     Departing  energy  is   the 
enthalpy  of  (1+    V^/^        )lb     of  products  and  that  part  of  the  chemical 
energy   (1-  7{coMb  )[C*/a.  ^jjJL         that  is  not  converted  into  molecular 

form.     Combustor  efficiency  as  used  here  is  defined  as: 


% 


COMl 


actual  energy  taken  on  by  the  gases 

energy  released  during  complete  combustion 
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The  resulting  equation  is 

hi  +  ^A^i  +  h^  -(\+^)/L4+(»-^B)Vv-f«  <2-12> 

Then 


ty*= 


kt-k; 


(2.13) 


The  lower  heating  value  of  the  fuel  q^   is  used  in  gas  turbine  analysis 
as  it  is  impractical  to  expand  the  gas  to  a  point  where  water  vapor  from 
the  products  of  combustion  will  condense  arid  contribute  any  energy  to 
the  cycle  (3). 

The  enthalpy  of  liquid  octane  with  respect  to  a  datum  of  octane 
vapor  having  h-0  at  absolute  zero  is  given  by 

^  -  .5V  1  *7  B  rU/^  Fu£L  (2. 14) 

The  enthalpy  of  other  hydrocarbon  combustibles  can  safely  be  taken  as 
the  same  (the  amount  of  energy  is  small)  (1). 

2.4  Turbine. 

If  the  turbine  is  considered  a  steady  flow  machine  with  kinetic 
energy  change  equal  to  zero,  the  works  of  the  fluid  for  adiabatic  flow 
are 

ideal:  W^-  (\*  V~W)(/ufx  "An)  <2'15> 

actual:  Wt   *(l  +  ^A)(l4  _  ^  (2.16) 

As  for  the  compressor,  the  foregoing  equations  can  be  obtained  from  the 
steady  flow  energy  equation.   The  turbine  efficiency  is  then 

7/  -  ACTUAL   FLUlQ  WQtiK      __       kgs  "  kr  (2   1?) 

<f"  IDEAL    WOAK  /L4^-^^ 
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2.5  Overall  Thermal  Efficiency. 

The  overall  thermal  efficiency  is  defined  as  net  work  out  divided 
by  energy  chargeable  against  the  cycle. 

-a        _We_ 

(fA-~~  Sc  (2.18) 

The  energy  chargeable  against  the  cycle  is  the  energy  released  by  the 
complete  combustion  of  the  fuel. 

The  equation  for  thermal  efficiency  then  follows  from  equations  (2.10), 
(2.16),  (2.18),  and  (2.19) 

^  ,  .  (l-vVV*Y^4~Ax)-  \kx~^\)  (2.20) 


2.6  Horsepower  Output. 

The  equation  for  horsepower  output  follows  directly  from  equations 
(2.18)  and  (2.19) 


hp-w^MW(-^Au4 


(2.21) 


where  m  has  units  of  lb  air  per  hour  and  2544  BTU  =  one  horsepower 
a  m 

hour. 


2.7  Specific  Fuel  Consumption. 

Specific  fuel  consumption  is  defined  as  lb  fuel  burned  per  unit 

n 

work  output.   If  the  horsepower-hour  is  used  as  the  unit  of  work  then 

specific  fuel  consumption  has  the  units  of  lb   fuel  per  horsepower-hour, 

m 


It  follows  that 


sfc  =^r —  (2-22) 
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3.   Regenerator  Surface  and  Core  Geometry. 

Certain  geometrical  relations  are  necessary  in  the  application  of 
the  basic  heat  transfer  and  flow-friction  data  to  a  design  problem.   One 
particular  set  and  form  of  these  relations  that  has  proven  convenient 
for  use  in  the  computer  program  is  given  here. 


MMMJ 


AAAAAA 


COLb* 

HOT 
COLD- 
HOT 
COib 


Figure  3  Simple  Plate -Fin  Counterflow  Regenerator  Core 

The  following  part  of  basic  data  are  given  for  all  surfaces  in  Kays 

and  London  (4) . 

b  -  plate  spacing 

<h  =  flow  passage  hydraulic  radius 

O  ~  ratio  of  total  transfer  area  of  one  side  of 

exchanger  to  volume  between  plates  of  that  side 

a  ■  plate  thickness 

In  addition,  the  ratio  of  fin  area  to  total  transfer  area  is  normally 

given. 

The  following  geometrical  factors  are  required  for  each  of  the  two 
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sides,  hot  and  cold,  of  the  complete  heat  exchanger  core. 

A  ■  total  transfer  area  of  one  side  of  exchanger 

A  ™  free  flow  area  of  one  side 
c 

Afr  ■  frontal  area  of  one  side 

L  =  flow  length  of  one  side 

V  =•  total  exchanger  core  volume 

°*  *  ratio  of  total  transfer  area  of  one  side  of  exchanger 
to  total  exchanger  volume 

CT»  ratio  of  free-flow  to  frontal  area  of  one  side  of 
exchanger 

V  ■  solid  core  volume 
s 

The  equations  below  give  the  relations  between  surface  and  core 
factors  for  one  side  of  the  exchanger.   Subscript  "1"  refers  to  any  one 
side  and  "2"  refers  to  the  other  side;  factors  without  a  subscript  are 
common  to  both  sides. 


(3.1) 


Ak  \  _  LA k),  - /. ^ \ 


(3.2) 


AA 


^^t±aJ£a _  (3.3) 


«,,A-,/    A      Wjr.  (3.4) 

1     V      1  L  A*.  J,   \  itk  /, 

^    AAi±tL.  (3#5) 


A.'-  (r/W\-  (AS.'te-l 
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A,-  (*LA#A 


(3.7) 


VmLtA^i-+A^-(Ae.+A4 


(3.8) 
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4.    Counterflow  Exchanger  Heat  Transfer  Theory 
4.1  Conduction  Heat  Transfer. 

When  a  temperature  gradient  exists  in  a  body  there  is  an  energy 
transfer  from  the  high  temperature  region  to  the  low  temperature  region. 
The  energy  is  transferred  by  conduction  and  the  heat  transfer  rate  per 
unit  area  is  proportional  to  the  normal  temperature  gradient. 


\"£ 


(4.1) 


When  the  proportionality  constant  is  inserted, 

where  "k"  is  called  the  thermal  conductivity  of  the  material, 


(4.2) 


k-L- 


T, 


(-»rX: 


dx 


A 


Figure  4  Elemental  Volume  for  One -Dimensional  Heat 
Conduction  Analysis 


If  the  steady  state  one-dimensional  case  is  considered 
which  integrates  to  give 

*  L 

providing  thermal  conductivity  is  constant. 


(4.3) 
(4.4) 


4.2  Convection  Heat  Transfer 


Figure  5  Sketch  Illustrating  Convection  from  a  Plate 
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The  convective  heat  transfer  rate  is  related  to  the  overall  temperature 
difference  between  the  wall  and  fluid  and  the  surface  area  as  follows: 


f=AA("fc-Tj) 


(4.5) 


where  "h"  is  called  the  convection  heat  transfer  coefficient.   It  may  be 
anticipated  that  convection  heat  transfer  will  have  a  dependence  on  the 
thermophysical  properties  of  the  fluid  (thermal  conductivity,  specific 
heat,  density,  viscosity),  the  speed  at  which  the  fluid  is  moving  and 
some  characteristic  dimension,,   If  a  dimensional  analysis  is  performed 
using  hi  u,  c  ,  {O   ,  k  ,  M.      and  L  as  variables,  the  variables  may  be 
arranged  into  three  independent  dimensionless  parameters.   The  para- 
meters depend  on  how  the  analysis  was  performed.   One  very  useful  set  of 
parameters  is 


k      f     P  U   l~        -D  .M  C<*  P 

^^TcT^'J  XT'  "faj  T—^W 


Pucf     wtJ 


(4.6) 


A  function  of  the  form 


t(&jfyRt)  =  o 


(4.7) 


must  therefore  exist.   Such  a  function  may  be  determined  experimentally 
in  the  form 

for  various  heat  transfer  surfaces  and  geometries.   The  generalized 
grouping  ^t   lr     is  defined  as  the  Colburn  j-factor  and  the  char- 
acteristic dimension  "L"  in  "Re"  is  designated  as  the  hydraulic  dia- 
meter, "D,  ". 
n 
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R*  = 


PlA   Pb   _.  &    D*  (4.9) 

A  large  number  of  curves  expressing  the  relationship  given  by  equation 
(4.8)  are  published  in  Kays  and  London  (1).   One  of  the  more  common  ex- 
perimental methods  for  obtaining  these  characteristics  is  described  in 
a  paper  by  Pucci,  Howard  and  Piersall  (4). 

After  determining  the  Colburn  j  factor,  St   is  easily  found  pro- 
viding the  fluid  properties  necessary  to  determine  Pr  are  obtainable. 


&=*M 


yi 


(4.10) 


The  convection  heat  transfer  coefficient  is  then  evaluated  by  using 


k^  vk  £(> 


(4.11) 


4.3  Fin  Effectiveness. 


t. 


Figure  6  Sketch  Illustrating  One -Dimensional  Conduction 
Through  a  Rectangular  Fin  with  Convection  at 
the  Surface 

Consider  the  one  dimensional  fin  exposed  to  a  surrounding  fluid  at 

a  temperature  Tqo   •   The  temperature  of  the  base  of  the  fin  is  T  . 

o 

Making  an  energy  balance  on  an  element  of  the  fin  of  thickness  dx  as 

shown  in  figure: 

(energy  in  left  face)=(energy  out  right  face)+(energy  lost  by  convection) 
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energy  in  left  face  =  #x  =  —  A  n~Tx~ 


energy  out  right  face  -  Jx+-U  .  _yA  ^  j^  =  -Jfi  (*J  +  ^J.^ 

Let  P  equal  the  perimeter  of  the  fin  then: 

energy  lost  by  convection  ■  <^4C  —  A  I  ^A  (  |  ""  loo) 
Combining  quantities,  the  energy  balance  yields 

7?-  if  (T-T-)-O  (*•!«) 


Let  (9  =  T  -  T^   .   Then  (4.12)  becomes 
One  boundary  condition  is 

a-  a,- To- To,  o***° 

For  the  case  of  a  fin  with  an  insulated  end,   dT/dx  =  0  and  d  0   /dx  = 

0  at  x  ■  Z. 

Let 

"*  -  ^-  (4.14) 

The  general  solution  for  (4.13)  may  be  written 

Solving  for  constants  C.  and  C  the  solution  is  obtained  as 

-#-  =  * +  i (^l6> 

The  total  heat  lost  by  the  fin  must  equal  the  energy  conducted  through 
the  base  of  the  fin 
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Then 


Fin  efficiency  is  defined  as 


r[$  -   actual  heat  transferred 


heat  which  would  be  transferred  if  entire  fin  area  were  at  base 

temperature 


which  becomes 


^    V/PiA  Qo£U(4Hi)_  sfc+Ji  (***)  (4.i8) 

To  apply  equation  (4.18)  in  the  computer  program  it  is  necessary  to 
put  it  in  terms  of  the  heat  exchanger  core  geometry.   Referring  to  Figure 
6, 

A  *  f  L 

so  that 


It  is  seen  from  Figure  3  that  the  fins  of  the  heat  exchanger  are 
symmetrical  with  respect  to  a  plane  midway  between  the  plates  separat- 
ing the  fluids.   Therefore,  the  temperature  gradient  in  the  fin  at  this 
plane  is  zero  and  the  above  insulated  end  solution  applies. 

Noting  that  Z  ?&    b/2,   the  fin  efficiency  becomes 

4.4  Surface  Effectiveness. 

The  total  effectiveness  for  a  finned  surface  may  be  defined  as  the 
ratio  of  the  total  heat  transfer  of  the  combined  area  of  the  surface  and 
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fins  to  the  heat  which  would  be  transferred  if  this  total  area  were 


maintained  at  the  base  temperature  T  .   Let  "A"  equal  the  total  heat 
transfer  surface  including  the  fin  area  Af  then 


^.•l-  ^MA 


(4.20) 


4.5  Overall  Coefficient  of  Heat  Transfer. 


T« 


T  T~ 

-WV — J_ /vW — £ — 'WV 


T 


^A 


X-A 


Figure  7  Sketch  Illustrating-One  Dimensional  Conduction 
and  Convection  through  a  Plane  Wall 

Consider  the  plane  wall  shown  above  exposed  to  a  hot  fluid  on  one 
side  and  a  cooler  fluid  on  the  other  side.   The  steady  state  heat  trans- 
fer is  expressed  by 


f=X»A(r»-i)-AA 


A* 


tt-T^-JAfa-Z) 


(4.21) 


The  heat  transfer  process  may  be  represented  by  the  resistance  network 
shown  and  the  overall  heat  transfer  is  calculated  as  the  ratio  of  the 
overall  temperature  difference  to  the  sum  of  the  thermal  resistances. 

Th-Tc 


^a+h-4-  y^A 


(4.22) 


The  overall  heat  transfer  coefficient  "U",  is  defined  by  the  re- 
lation 

4~U  A  ATovz <all  (4.23) 
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where  A  is  some  suitable  area  for  heat  flow.   Then 


y*»  +  ¥  +  '//< 


A  fT>l 

Now  consider  a  plate  of  thickness  "a"  with  extended  surfaces  or  fin 
areas  on  either  side.   The  overall  heat  transfer  coefficient  becomes 


U- 


-r 


a 


(4.24) 


-V- 


V  -*  H         (A"/A«  )  -4      '  (Ac/A  M  )  %,C  /c 
where  U  is  based  on  a  unit  of  hot  side  total  area  including  fin  or  ex- 
tended area.  A  denotes  the  average  wall  area.   In  the  computer  program 
w 

the  very  small  wall  resistance  is  neglected  in  calculations  so  that 


u 


v^ 


+ 


(kM*l\tAt 


4.6 


Exchanger  Heat  Transfer  Effectiveness 
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FLOW  L.  E/V6TH 

Figure  8  Temperature  Profiles  in  Counterflow 
Heat  Exchanger 

The  actual  heat  transfer  in  a  heat  exchanger  may  be  computed  by 

calculating  either  the  energy  lost  by  the  hot  fluid  or  the  energy  gained 

by  the  cold  fluid.   Referring  to  Figure  8, 


Heat  exchanger  effectiveness  is  defined  as 


£   = 


actual  heat  transfer 


maximum  possible  heat  transfer 
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To  determine  the  maximum  possible  heat  transfer  for  the  exchanger  it  is 
necessary  to  recognize  that  this  maximum  value  could  be  attained  if  one 
of  the  fluids  were  to  undergo  a  temperature  change  equal  to  the  maximum 
temperature  difference  present  in  the  exchanger,  (  ~\^/fj—    Tc/ ///    )• 
The  fluid  which  might  undergo  this  maximum  temperature  difference  is  the 
one  having  the  minimum  value  of  /yk-C-jo      since  the  energy  balance  re- 
quires that  the  energy  received  by  one  fluid  is  equal  to  that  given  up 
by  the  other  fluid.   If  the  fluid  with  the  larger  value  of  rt*>£fO     goes 
through  the  maximum  temperature  difference,  the  other  fluid  would  under- 
go a  temperature  difference  greater  than  maximum.   The  maximum  heat 
transfer  is  therefore  expressed  as 


f»**  =  (*£,)»,„   (TTyv-Tc,^)  (4.26) 


The  "minimum"  fluid  may  be  either  the  hot  or  cold  fluid  depending  on  the 
mass  flow  rates  and  specific  heats. 

£  _  ( ^ cp) h  (Thx - TV>) o ut  I   _  [a^£/>)c  (Ty^r  -Ti,/*/ )       (4.27) 

An  expression  for  effectiveness  in  counterflow  will  now  be  derived.   Re- 
ferring to  Figure  8,  the  heat  transferred  through  an  element  of  area  dA 
may  be  written 

The  heat  transfer  could  also  be  expressed  as 

*%   =U(Th-ToVA  (4.29) 

Equations  (4.28)  and  (4.29)  combine  to  give 
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Th  -  Tc 


Integration  yields 


It  \ 

{   'H,60T-  Tc, /V  J 


c 


c 


or 


fT. 


UA 


H>o  yT 


If    £vi<  Cc      or  C»*C»,4 


/W 


Then 


and 


T"  ^ 


From   (4.32) 


(4.30) 


(4.31) 


(4.32) 


(4.33) 


(4.34) 


(4.35) 


-  £ 


-  £ 


t 


_c£\\J£ 

Cc  J  ^* 


(4.36) 
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Solving  for   £ 


€-■ 


If    Lc^G  or   Cc  -  L^ 


l  A>       , 


the   same  procedure  yields 


In  either  case 


U      c«  J  Cc 


€  = 


C*>d 


iM 


C/nA*J    Ca>'A/ 


£  = 


-  £ 


\}        CaiaxJ' 


VTU 


^  m  / 


v     - 


C/waX  £ 


C*AXj 


(4.37) 


(4.38) 


(4.39) 


(4.40) 


where 


4/rU-^b 


is  defined  as   the  number  of  heat   transfer  units 
possible   to  express 


In  general   it  is 


£  =  p  ^/SlJ^LjU^L    ,  Flow  arrangement) 


C* 


AX 
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5.   Exchanger  Flow-friction  Theory. 


Figure  9  Heat  Exchanger  Core  Model  for  Pressure-Drop 
Analysis 

5.1  Core  Friction. 

Friction  factor  is  defined  on  the  basis  of  an  equivalent  shear 
force  in  the  flow  direction  per  unit  of  wetted  area.   For  most  of  the 
surfaces  considered  it  is  a  combination  of  viscous  shear  and  pressure 
force  but  no  attempt  is  made  to  separate  these  two  effects.   The  pres- 
sure force  (form  drag)  is  most  significant  in  interrupted  or  wavy  fin 
surfaces  where  fluid  separation  from  these  surfaces  occurs.  Thus 


or 


(5.1) 


(5.2) 


It  follows  directly  then  that 


3  Ac    *-f   yAc  13^ 


(5.3) 


where   /^   equals  the  average  density  in  the  core  and  "G"  is  based 
on  free  flow  area.   Curves  showing  friction  factor  as  a  function  of 
Reynolds  number  are  published  in  Kays  and  London  (1).   They  are  plotted 
for  the  same  surfaces  as  the  Colburn  j -factor  mentioned  in  the  previous 
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section. 

5.2  Flow  Acceleration. 

A  term  accounting  for  flow  acceleration  must  also  be  included  be- 
cause of  the  change  in  density  as  the  fluid  flows  through  the  core.  A 
force  balance  using  change  in  momentum  yields 

F»  /^(?<i-- W)  (5.4) 


or 

The  core  friction  term  and  the  flow  acceleration  term  may  now  be  com- 
bined to  give  the  pressure  drop  between  points  "a"  and  Mb". 

ip^fA^+^i/iA  (56) 


5.3  Entrance  and  Exit  Effects. 

A  typical  gas  turbine  regenerator  installation  usually  involves  a 
flow  contraction  at  the  heat  exchanger  core  entrance  and  a  flow  expan- 
sion at  the  core  exit.   They  introduce  additional  flow  stream  pressure 
drops.   This  problem  has  been  investigated  both  analytically  and  experi- 
mentally by  Kays  (5) . 

The  entrance  pressure  drop  is  made  up  of  two  parts.   The  first, 
(P-  -  P  ).,  is  the  pressure  drop  which  would  occur  due  to  a  flow  area 
change  alone.   The  second,  (P  -  P  )  ,  is  the  pressure  loss  due  to  the 
irreversible  free  expansion  that  always  follows  the  abrupt  contraction. 

First  considering  the  pressure  drop  due  to  area  change  alone, 

fi+  if  «|l-a  +  -tf  VLir  (5.7) 
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The  pressure  change  from  "1"  to  Ma"  is  very  small  relative  to  the  total 
pressure  so  that  p  &   z?  .   Therefore 

AfrU^AcU*  (5.8) 

Substitution  of  equation  (5.8)  into  (5.7)  yields 

The  analysis,  as  presented  in  Kays,  of  the  pressure  drop  due  to 

free  expansion  is  somewhat  more  complex  and  only  a  brief  description  of 

the  equations  utilized  in  the  computer  program  will  be  given  here.   The 

contraction  loss  coefficient,  K  is  defined  as  follows: 

c 


(5.10) 


K  can  be  evaluated  from 
c 


^S  l-J-^ciTL.  ^^  I,  (5#u) 


K,  is  a  function  of  Reynolds  number  and  flow  geometry.   For  laminar 

flow 

K,  =  1.333  (circular  tubes)  (5.12) 

K,  ■  1.390  (square  tubes) 

K,  ■  1.430  (triangular  tubes) 

For  turbulent  flow  in  circular  tubes 

Ak»/.o*0£ff(4£,)  +  .ormif^T-H  (5.13) 


where 


and 


f-     0.0W  Re  (5.14) 


k/(swMe}-\        -J  iq  ICi((TR\Att*Lr}-l       _(  ^   (5.15) 
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C  is  a  function  of  the  contraction  ratio  0~  and  is  presented  graphical- 
ly for  circular  tubes  and  a  one -dimensional  contraction.   Kays  recommends 
that  the  one-dimensional  contraction  be  used  in  evaluating  triangular 
and  square  tubes.   The  curves  of  C  vs.  cj"  were  fit  with  straight  line 
approximations  for  use  in  the  computer  program  and  are  shown  in  Ap- 
pendix B. 

The  exit  pressure  rise  is  similarly  made  up  of  two  parts.   The 

first,  (P,  -  P„),,  is  the  pressure  rise  which  would  occur  due  to  area 
b    l   i 

change  alone.  The  second,  (P  -  P9)9»  is  the  pressure  loss  associated 
with  the  irreversible  free  expansion  and  momentum  changes  following  an 
abrupt  expansion.   The  first  term  follows  from  the  same  analysis  as  the 
corresponding  term  in  the  entrance  pressure  drop. 

The  free  expansion  term  is  given  by 

■X- 

(Pi-R-k  „  If    %C_  ,.  U    _6_  (5.17) 

where 

to-  \-ZKc(<r+  r^  (5.i8) 

and  K,  is  evaluated  using  (5.12),  (5.13)  and  (5.15). 
a 

K  and  K  are  presented  graphically  in  Kays  and  London  for  the 
c      e 

circular,  triangular  and  square  tubes j  however,  the  preceding  equations 
were  used  for  evaluation  in  the  computer  program. 

If  the  core  friction,  flow  acceleration,  and  the  entrance  and  exit 
effects  are  now  combined 
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(?-R-4J[6-'^M^-'V*^l  +  (cr^ic)^ 


3-t 


O^VfeVl^^^fcf+^l  +  Kt^l  C»-l» 
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6.    Computer  Program  Development. 
6.1  General  Approach. 

As  discussed  in  the  introduction,  the  program  simulates  a  gas  tur- 
bine arrangement  running  under  constant  load  conditions.   The  nature  of 
the  problem  necessitates  that  an  iterative  approach  be  taken  in  solving 
for  most  of  the  desired  outputs.   The  program  is  flexible  in  that  six 
different  sets  of  output  data  are  obtainable  depending  on  what  type  of 
information  the  user  desires. 

Program  inputs  can  be  divided  into  four  categories  as  follows: 

a.  known  information  describing  the  cycle  as  a  whole 

b.  information  describing  regenerator  core  characteristics 

c.  assumed  quantities  necessary  to  start  the  solution 

d.  special  inputs  depending  on  the  type  of  output  desired 
The  first  three  categories  are  referred  to  hereafter  as  general  inputs. 

Program  outputs  can  be  divided  into  two  categories  as  follows: 

a.  special  outputs  depending  on  which  set  of  output  data  is 
desired 

b.  general  outputs  common  to  all  sets  of  output  data 

The  program  is  executed  largely  through  the  use  of  subroutines  so 
that  characteristics  of  various  gas  turbine  arrangements  are  easily 
interchangeable.  All  general  outputs  are  solved  for  in  one  rather 
long  subroutine  called  XMAIN.   Subroutine  XMAIN  utilizes  all  the  other 
subroutines  as  the  solution  progresses.   Subroutine  XMAIN  is  described 
in  detail  later  in  this  section;  however,  the  general  approach  for  the 
solution  of  the  various  unknowns  in  the  gas  turbine  cycle  is  as  fol- 
lows: 

a.   Referring  to  Figure  1,  the  temperature  and  pressure  at 
point  "2"  are  solved  for  knowing  the  temperature  and  pressure  at  point 
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"1",  the  compression  ratio  and  the  compressor  efficiency. 

b.  The  pressures  at  points  "5"  and  "3"  are  computed  using 
assumed  values  for  pressure  drops  through  the  regenerator. 

c.  An  estimate  of  the  temperature  at  point  "5"  is  computed 
using  ^-t   and  the  compression  ratio. 

d.  The  temperatures  at  points  "3M  and  "6"  are  computed  using 
the  estimated  temperature  at  point  "5"  and  an  assumed  value  for  re- 
generator effectiveness. 

e.  Better  values  are  computed  for  the  regenerator  pressure 
drops  using  theory  from  section  5  and  an  assumed  mass  flow  rate. 

f.  Better  values  are  computed  for  the  temperatures  at  points 
"5",  "3"  and  "6". 

g.  A  new  value  is  computed  for  regenerator  effectiveness 
using  the  theory  from  section  4  and  fed  back  into  step  "d". 

h.   The  process  is  repeated  until  all  temperatures  and  pres- 
sures converge. 

i.   A  new  value  of  mass  flow  rate  is  now  computed  by  comparing 
the  calculated  value  of  power  output  to  the  desired  power  output. 

j.   The  entire  process  is  then  repeated  until  mass  flow  rate 
converges. 

k.  All  remaining  general  outputs  are  now  computed. 

6.2  General  Inputs. 

a.  The  following  inputs  describe  the  cycle  as  a  whole: 
T.        compressor  inlet  temperature 
P-        compressor  inlet  pressure 


P_  /P     compressor  pressure  ratio 
'(c.  compressor  efficiency 
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T,  turbine  inlet  temperature 

p  hot  side  regenerator  exhaust  pressure 
6 

\-t  turbine  efficiency 

'(comb.  combustor  efficiency 

q relative  pressure  drop  through  combustor 

HP  net  horsepower  output 

Tf  fuel  temperature  at  combustor  inlet 

q^  lower  heating  value  of  fuel 

R  gas  constant  for  air 
b.   The  following  inputs  describe  the  regenerator  core  character- 
istics: 

/O  density  of  regenerator  material 

$  cold  side  heat  transfer  area/volume  between 
cold  side  plates 

ft  hot  side  heat  transfer  area/volume  between 
hot  side  plates 

J^i  cold  side  plate  spacing 

jlr^.  hot  side  plate  spacing 

Y\c  cold  side  hydraulic  radius 

«*-A  hot  side  hydraulic  radius 

oc  cold  side  fin  thickness 


u 


hot  side  fin  thickness 


A_,.      cold  side  fin  or  extended  area/total  cold  side 
transfer  area 

Aj  t      hot  side  fin  or  extended  area/total  hot  side 
transfer  area 

a        plate  thickness  between  cold  and  hot  sides 

k        fin  thermal  conductivity 

The  following  are  assumed  input  quantities: 

m        mass  flow  rate  of  air 
a 
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3 Atr  relative  pressure  drop  through  cold  side  of  re- 

y\_  generator 

(, Af  P  relative  pressure  drop  through  hot  side  of  re- 

Pp  generator 

£  regenerator  effectiveness 


6.3   Subroutines  Other  Than  Subroutine  XMAIN 

a.  Subroutine  KEKC  (SIGMA,  RE,  XKE,  XKC) 

Input  arguments  are   (J~   and  Re  ;  computes  Ke  and  Kc 
using  equations  from  section  5. 

b.  Function  PRAND  (T) 

Input  argument  is  T  ;  computes  Prandtl  number  of  air  from 
equation  obtained  by  fitting  curve  to  data  given  in  Hilsenrath  (6). 
See  Appendix  B  for  data. 

c.  Function  VIS  (T) 

Input  argument  is  T  ;  computes  /A       of  air  from  equation 
obtained  by  fitting  curve  to  data  given  in  Hilsenrath  (6).   See  Ap- 
pendix B  for  data. 

d.  Function  SPEC  (A) 

Input  argument  is  T   ;  computes  c  of  air.   The  equations 
used  for  computation  in  this  and  the  next  four  subprograms  are  given  in 
Howard  (2). 

e.  Function  TEMR  (B) 

Input  argument  is  Pr  ;  computes  T  of  air.   Pr  as  used  here 
refers  to  relative  pressure,  not  to  be  confused  with  Prandtl  number. 

f.  Function  ENTH  (B) 

Input  argument   is     T    ;    computes  h  of  air. 

g.  Function  TEMH  (B) 

Input  argument  is  h;  computes  T  of  air. 
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h.   Function  PREP  (B) 

Input  argument  is  T;  computes  Pr  of  air. 

i.   Function  FUENTH  (T) 

Input  argument  is  Tf;  computes  enthalpy  of  fuel  from  equation 

given  in  section  2. 

j.    Function  FRICC  (RE) 

Input  argument  is  regenerator  cold  side  Re;  computes  cold 

side  f  from  equation  obtained  by  fitting  curve  to  experimental  data. 

See  Appendix  B  for  details  on  f  and  j  vs.  Re  data  used  in  this  and  the 

following  three  subprograms. 

k.   Function  FRICH  (RE) 

Input  argument  is  regenerator  hot  side  Re;  computes  hot  side 

f  from  equation  obtained  by  fitting  curve  to  experimental  data. 

1.   Function  STPRC  (RE) 

Input  argument  is  regenerator  cold  side  Re;  computes  cold  side 

j  from  equation  obtained  by  fitting  curve  to  experimental  data. 

m.   Function  STPRH  (RE) 

Input  argument  is  regenerator  hot  side  Re;  computes  hot  side 

j  from  equation  obtained  by  fitting  curve  to  experimental  data. 

n.   Function  GENFX  (TU,  RATIO) 

Input  arguments  are  C  .  /C    and  NTU;  computes  £-       from 

mm  max 

equation  developed  in  section  4. 
o.    Subroutine  XWRITE 

This  subroutine  contains  all  general  output  print  statements, 
p.   Subroutine  XOPEN 

Computes  open  cycle  outputs 
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6.4  General  Outputs. 

t> 

m  mass  flow  rate  of  air 
a 

€.  regenerator  effectiveness 

Vf/a.  fuel  to  air  ratio 

SFC  specific  fuel  consumption 

V  regenerator  core  volume 

Wt.  regenerator  core  weight 

P_  compressor  discharge  pressure 

T  compressor  discharge  temperature 

P  cold  side  regenerator  outlet  pressure 

T,  cold  side  regenerator  outlet  temperature 

P,  turbine  inlet  pressure 

P  turbine  exhaust  pressure 

T,.  turbine  exhaust  temperature 

T,  regenerator  hot  side  outlet  temperature 
o 

Re  regenerator  cold  side  Reynolds  number 
c 

Re  regenerator  hot  side  Reynolds  number 

NTU  number  of  exchanger  heat  transfer  units 
Cmin/Cmax  capacity-rate  ratio 

/aP\  -  3A*P  .  4A3P  ,  «ArP 


6.5  Detailed  Developmentof  Subroutine  XMAIN 

Assumptions  made  in  the  development  of  subroutine  XMAIN  are: 

a.   Combustion  products  have  a  negligible  effect  on  the 
fluid  properties  of  gases  leaving  the  combustor. 

b-  f[d      ,  '[ ±    ,  '?Ct)M§       »  and  — pr remain  constant  dur- 

ing  the  iteration  process. 
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and 


The  first  step  in  the  subroutine  is  to  calculate  T?  from  T. ,  P. 


P.   ' 


k^tiiTi)  (6.1) 


R  , 

hiff=V(TW)  (6.5) 

JL=  <^(A^)  (6.7) 


Next,  h.  and  P   .are  computed  from  T. . 
*  4      rel,      r         4 


h^0(\)  (6.8) 


3AZ  P 
P  and  P  are  now  computed  using  the  assumed  values  for  -=j=j and 

*ArP 

P'  K«f?(^  <6-10> 


(3«R-(i^E)p,  (6.11) 

S-R/O'-^)  (6-12) 


It  next  becomes  necessary  to  get  an  estimate  of  T_.   This  is  done  by  using 
a  cycle  with  no  regenerator  and  no  pressure  drop  through  the  combustor. 
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(6.13) 


fc^-R*  (*■)-&«•(£)        (6-i4) 

~Trs  =  ^(R-^s)  (6-15) 

L's=^(T^  (6.16) 

l^-^(Kr)  (6.i8) 

T„  and  T,   are  now  computed  using  the  value  just  computed  for  T_  and 
3      o  D 

and  assumed  value  of  £ 

Y3-  fe(rr-Tt.)-v-Ti.  (6-19) 

T*  -  Tr -^  (Tj—T^)  <6-20> 

Temperatures  and  pressures  at  points  "2",  "3",  "5"  and  "6"  are  now 
known  or  have  been  estimated.   Better  values  are  computed  for  the  pres- 
sure drops  through  the  regenerator  using  the  theory  from  section  5,  an 
assumed  value  of  m  and  fluid  properties  computed  from  the  temperatures 
and  pressures  at  the  above  four  points.   For  the  regenerator  cold  side: 

(6.21) 

(6.22) 

(6.23) 

(6.24) 


Vs 

■R(T,/RJ 

v2  = 

•  R  (Vs) 

T 

.  T.+  Tj 

Z 

V~>t 

2- 
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V-Lx  A«v 


(6.25) 


°<c- 


(6.26) 

(6.27) 
(6.28) 
(6.29) 
(6.30) 
(6.31) 
(6.32) 
(6.33) 


p-P,  6-c 


(i-aeVlci,0^(lfe-«)  +  fe^^ 


Ac/     * 


*-(<s*-H-K*)&] 


(6.34) 


The  pressure  drop  through  the  regenerator  hot  side  is  computed  in  the 
same  manner  except  that  the  fuel  to  air  ratio  must  be  included  in  cal- 


culating G,. 


where 


(6.35) 


ty*= 


(k*-kz) 


kir+>(a>«bp-k4) 


(6.36) 
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P»  and  P  are  now  recalculated  using  the  results  just  obtained. 


B-R-(S-fO  (6.37) 

R-B+(^-R)  (6-38) 


A  new  value  of  T,.  is  computed  as  follows: 


R-R-(4^P)P3  (6.39) 

iwirfB  P^^(^)  (6-40) 

"!}**  ^(f?*^  <6-41) 

krf*^(T?0  (6.42) 

Kf-  /m-'^/h-ZmV)  (6.43) 

TV  -  ^  (A.J-)  (6.44) 

Using  this  value  of  Tc  and  the  assumed  £   ,  new  values  for  T0  and  T, 

j  Jo 

are  computed.      T^   \^y  Mt    .   Mk      >  h3>       f  */*.     >   V   ReC  and  Reh 

are  then  recomputed  using  the   above  values   of  T0   and  T-. 

J      o 

A  better  £   can  now  be  determined  using  the  theory  from  section 
4. 

fa*  </>(R*j*)  (6.45) 

jt~   ^(R«,c)  (6.46) 

R>S^(W)  (6-47) 

R-/*^(T^^  (6.48) 

S^A-    ^k/?v//d  (6-49) 
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£fc.=  v^<c  (6'50) 

tfl,*-   ^  (T^,k)  (6.51) 

C/p/C  =  (^  (T^c)  (6.52) 

^s    6/>(Sa<:^A  (6.53) 

^  *S^c6tCAc  (6.54) 


C 


(6.55) 


(6.56) 


^A=  K4e*(i-ft,*)  (6.5?) 

^/C-   1-A*>c(l-^C)  (6.58) 

l> j (6.59) 

LA.  =    ^*(l+lfy*)<V,A.  (6.60) 


c  =  >*•*•*»(.«,£  (6.61) 

(6.62) 


(6.63) 
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The  above  value  of  €.       is  now  fed  back  into  equation  (6.19)  and  the 
entire  set  of  calculations  from  (6.19)  through  (6.63)  is  repeated  until 
the  desired  degree  of  convergence  is  obtained.  Many  different  criteria 
for  determining  convergence  could  be  used;  however,  numerous  trials  in- 
dicate that  the  regenerator  cold  side  pressure  drop  converges  at  a  slower 
rate  than  any  of  the  other  computed  outputs.  When  the  value  of  (P  -  P-)  , 
computed  from  equation  (6.34),  converges  to  within  one  percent  of  its 
previous  value  the  subroutine  leaves  the  above  loop  and  proceeds  to  the 
next  step. 

Thermal  efficiency  and  the  true  value  of  mass  flow  rate  needed  to 
give  the  desired  horsepower  output  at  that  thermal  efficiency  are  now 
computed. 


Then 


2-^44 


(6.65) 


where  rn^is  the  assumed  value.   The  above  calculated  value  of  horsepower 
is  then  compared  with  the  desired  value.   If  the  two  agree  within  .1 
percent  of  the  desired  horsepower  the  subroutine  ends;  if  they  do  not  a 


new  value  of  m  is  calculated 
a 


,  .  2J-4  4  HP 
/*"«. — „ — (6-66> 

ntk  If}/*  fy 


where  HP  is  the  desired  horsepower  output.   This  value  of  m  is  then 
fed  back  into  the  subroutine  where  the  process  of  calculating  pressure 
drops,  effectiveness,  thermal  efficiency  and  horsepower  is  repeated 
until  horsepower  converges.   Convergence  is  normally  rapid  and  two  or 
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at  most  three  times  through  the  set  of  calculations  is  sufficient. 

The  subroutine  then  proceeds  to  calculate  the  remainder  of  the 
general  outputs. 

SFC  ^  Zr41  (6.67) 

Vc=   L[Afv-(Ac^+Acye\]  <6-6*> 

Wt.   «    /^X\/<r  (6.69) 

/^=^P  +  ^P^^£.  (6.70) 

V  r   /total    ^  ft  ™ 

6.6  Special  Inputs  and  Outputs. 

As  mentioned  before,  the  program  is  flexible  in  that  different  sets 
of  output  data  are  obtainable  depending  on  what  type  of  information  the 
user  desires.   The  five  possible  sets  of  output  and  corresponding  inputs 
are  described  below. 

a.  Output  type  one 

In  some  cases  the  frontal  area  and  regenerator  core 
length  may  be  fixed  at  some  predetermined  values.   The  user  may  then 
determine  gas  turbine  performance  for  that  frontal  area 

and  length  using  output  type  one.   The  frontal  area  and  length  along  with 
the  general  inputs  are  read  into  subroutine  XMAIN  and  all  general  out- 
puts are  computed. 

Inputs:   general,  length,  frontal  area 

Outputs:  general 

b.  Output  type  two 

It  may  be  desirable  to  fix  the  overall  thermal  efficiency 
at  a  predetermined  value.   There  may  also  be  some  latitude  possible  in 
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choosing  the  frontal  area  of  the  regenerator  providing  the  length  re- 
mains between  specified  limits.   Output  type  two  is  used  to  determine 
a  frontal  area  for  each  length  under  consideration  that  will  give  the 
desired  overall  thermal  efficiency.   The  shortest  length  under  considera- 
tion, the  longest  length  under  consideration,  the  size  of  the  increments 
desired  between  the  shortest  and  longest  lengths,  the  desired  thermal  ef- 
ficiency and  all  general  inputs  are  read  into  the  computer  program.   Sub- 
routine XMAIN  is  then  called  with  the  shortest  length  and  an  initial 
value  of  frontal  area,  determined  by  the  computer,  as  arguments.   Frontal 
area  is  then  incremented  with  length  held  constant  until  the  computed 
thermal  efficiency  is  in  agreement  with  the  desired  thermal  efficiency. 
The  same  procedure  is  repeated  for  all  other  lengths  under  considera- 
tion. 

Inputs:   general,  shortest  length  under  consideration, 
longest  length  under  consideration,  size  of  increments  desired  between 
shortest  and  longest  lengths,  desired  thermal  efficiency  . 

Outputs:  general,  required  frontal  area  for  each  length 
considered. 

c.   Output  type  three. 

Output  type  three  is  similar  to  two  in  that  the  overall 
thermal  efficiency  is  set  at  a  predetermined  value.   In  this  case,  how- 
ever, it  is  assumed  that  there  is  some  latitude  possible  in  choosing  the 
length  providing  the  frontal  area  remains  between  specified  limits.  A 
length  must  therefore  be  determined  for  each  frontal  area  under  considera- 
tion that  will  give  the  desired  thermal  efficiency.   The  same  computing 
procedure  is  followed  as  in  output  type  two  except  that  frontal  area  is 
held  constant  and  length  is  incremented  until  the  desired  thermal  effi- 
ciency is  obtained. 
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Inputs:   general,  smallest  frontal  area  under  considera- 
tion, largest  frontal  area  under  consideration,  size  of  increments  de- 
sired between  smallest  frontal  area  and  largest  frontal  area,  desired 
thermal  efficiency. 

Outputs:  general,  required  length  for  each  frontal  area 
considered. 

d.  Output  type  four. 

In  some  cases  it  may  be  desirable  to  find  the  minimum 
frontal  area  that  will  allow  a  predetermined  thermal  efficiency.   In 
other  words  make  the  frontal  area  as  small  as  possible  without  giving 
consideration  to  length.   Output  type  four  is  used  to  determine  this 
minimum  frontal  area  and  the  length  that  must  be  used  with  it  to  obtain 
the  desired  thermal  efficiency. 

All  general  inputs  along  with  the  desired  thermal  efficiency  are 
read  into  the  program.   Subrouting  XMAIN  is  called  and  an  initial  length, 
depending  on  the  mass  flow  rate  and  internal  core  geometry,  is  deter- 
mined by  the  computer.   Frontal  area  is  then  incremented  until  the  de- 
sired thermal  efficiency  is  obtained.  A  new  length  is  then  selected  and 
again  frontal  area  is  incremented  until  the  frontal  area  corresponding  to 
the  new  length  gives  the  desired  thermal  efficiency.   This  process  is 
repeated  and  frontal  areas  are  compared  until  a  minimum  allowable  front- 
al area  is  obtained. 

Inputs:  general,  desired  thermal  efficiency 
Outputs:  general,  minimum  frontal  area  for  desired 
thermal  efficiency,  length  to  be  used  with  that  frontal  area 

e.  Output  type  five. 

The  user  may  find  it  helpful  to  determine  the  maximum 
thermal  efficiency  possible  with  a  fixed  allowable  frontal  area  or  range 
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of  frontal  areas.   Output  type  five  computes  the  maximum  thermal  ef- 
ficiency and  the  length  that  must  be  used  with  the  frontal  area  under 
consideration  to  obtain  that  thermal  efficiency.  The  frontal  area  under 
consideration  is  read  into  subroutine  XMAIN  and  length  is  incremented 
until  thermal  efficiency  reaches  a  maximum. 

Inputs:   general,  smallest  frontal  area  under  considera- 
tion, largest  frontal  area  under  consideration,  size  of  increments  de- 
sired between  smallest  frontal  area  and  largest  frontal  area. 

Outputs:  general,  length  corresponding  to  each  frontal 
area  considered  for  maximum  thermal  efficiency. 

f.  Output  type  six. 

Output  type  six  is  utilized  in  determining  thermal  ef- 
ficiency, mass  flow  rate,  specific  fuel  consumption,  fuel  to  air  ratio, 
compressor  discharge  temperature  and  pressure  and  turbine  exhaust  tempera- 
ture of  the  simple  open  cycle.   This  information  can  be  used  in  compari- 
son with  other  output  types  in  setting  a  standard  for  regenerator  per- 
formance . 

Inputs:  general  except  regenerator  characteristics 

Outputs:  as  given  above 

g.  Tabular  summary  of  special  inputs  and  outputs 


Output 
Type 

fr 

fr,min 

fr,max 


One 

I 


Two 

0 
0 


Three 

I 

I 


Four 


Five 

I 
I 


Six 


mm 


max 


opt 

(tk,<hAX 
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GENERAL 
INPUTS 


SPECIAL 
INPUTS 


KE,  KC 


*.  J 


FLUID 
PROPERTIES 


REGENERATOR 
EFFECTIVENESS 


MAIN  PROGRAM 


SUBROUTINE 
XMAIN 


SPECIAL 
OUTPUTS 


SUBROUTINE 
XWRITE 


GENERAL 
OUTPUTS 


FIGURE  10.   FLOW  CHART  OF  COMPUTING  PROCEDURE 
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7.   Results. 

The  ORENDA  OT-4  gas  turbine  was  utilized  as  a  model  in  the  develop- 
ment and  testing  of  the  computer  program.   This  engine  meets  a  require- 
ment that  was  generated  by  the  U.  S.  Army  and  Navy  in  1960  and  is  dis- 
cussed in  a  paper  by  Quan  (7).   The  design  target  at  600  hp  on  specific 
fuel  consumption  was  .40  lbm  of  fuel  per  hp  hr,  which  is  equivalent  to 
a  thermal  efficiency  of  34.6%.  The  pressure  ratio  ranges  from  four  to 
five.   The  following  were  additional  design  targets: 

compressor  efficiency  86% 

turbine  inlet  temperature  2200 °R 

turbine  efficiency  89% 

combustor  efficiency  98% 

regenerator  effectiveness  82.5% 

A  few  additional  quantities  were  required  to  complete  the  general  inputs 
in  the  computer  program.  Their  values  were  selected  as  realistically  as 
possible  as  follows: 

compressor  inlet  temperature  560°R 

-2 
compressor  inlet  pressure  2120  lbf-ft 

fuel  inlet  temperature  560°R 

lower  heating  value  of  fuel  19,000  BTU-lbm" 

-2 
regenerator  hot  side  exhaust         2120  lbf-ft 

pressure 
Six  cores  were  selected  at  random  for  evaluation  and  incorporated 
into  the  above  gas  turbine  arrangement.   The  characteristics  of  the 
first  five  shown  on  the  following  pages  are  from  Kays  and  London  (1); 
the  last  one  from  Pucci  et.al.(8).  The  friction  factor  and  Colburn 
j -factor  data  utilized  in  obtaining  functional  relationships  for  each 
core  is  listed  in  Appendix  B  along  with  the  curve  fitting  results.   The 
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cores  were  made  of  stainless  steel,  aluminum  or  nickel;  however,  the 
density  of  stainless  steel  was  used  as  the  input  in  the  program  for  all 
cores  to  obtain  a  common  basis  for  weight  comparison.  Output  types  one 
through  six  were  computed  for  each  core.   The  desired  thermal  efficiency 
of  .34  was  set  for  output  types  two,  three  and  four.  Required  frontal 
areas  were  computed  to  match  lengths  of  3.0  ft  through  7.0  ft  in  output 

type  two  and  required  lengths  were  computed  to  match  frontal  areas  of 

2  2 

3.0  ft  through  7.0  ft  in  output  type  three.  All  core  incrementing  in 

output  type  three  was  done  in  .1  ft  increments  so  that  lengths  determined 

for  >\th.  ~   «34  are  only  accurate  to  within  .1  ft.  The  thermal  efficiency 

obtainable  with  these  computed  lengths  may  therefore  be  slightly  higher 

or  lower  than  .34.  The  characteristics  of  each  core  along  with  results 

selected  from  the  various  program  output  types  are  tabulated. 

The  tabular  results  are  followed  by  six  sets  of  curves  which  il- 
lustrate the  range  obtainable  in  overall  dimensions  and  weight  depend- 
ing on  the  core  selected.  The  set  of  curves  shown  in  Figure  12  illus- 
trates the  lengths  that  must  be  used  over  a  range  of  frontal  areas  to 
obtain  a  thermal  efficiency  of  .34.  Figure  13  gives  regenerator  core 
weights  as  a  function  of  frontal  area  to  obtain  a  thermal  efficiency  of 
.34  and  Figure  14  is  a  plot  of  regenerator  effectiveness  vs.  total  re- 
lative pressure  drop  for  a  thermal  efficiency  of  .34.  Figure  15  il- 
lustrates the  maximum  thermal  efficiency  obtainable  as  a  function  of 
frontal  area  and  the  optimum  lengths  for  each  frontal  area,  i.e.,  the 
lengths  necessary  to  obtain  maximum  thermal  efficiency  are  shown  in 
Figure  16.   Core  weights  as  a  function  of  frontal  area  for  maximum 
thermal  efficiency  are  shown  in  Figure  17. 

A  sample  of  the  six  computer  output  types,  generated  using  core  1, 
is  shown  following  the  curves.   Output  type  one  is  shown  on  page  73, 
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two  on  pages  74  through  76  ,  three  on  pages  77  through  79  ,  four  on 
page  80  ,  five  on  pages  81,  through  83  and  six  on  page  84  . 

As  mentioned  above,  the  six  cores  evaluated  were  selected  at  ran- 
dom.  No  systematic  attempt  was  made  in  finding  an  optimum  core  from 
the  standpoint  of  external  dimension  or  weight  considerations.   The 
results  serve  only  to  illustrate  the  usefulness  of  the  computer  program 
in  selecting  the  best  core  for  a  given  gas  turbine.  A  few  very  general 
observations  regarding  the  geometries  investigated  are  as  follows: 

1.  The  desired  thermal  efficiency  was  not  obtainable  with  the 
fine  mesh  cores  in  the  small  frontal  area  range,  regardless  of  length. 

2.  Significantly  longer  lengths  were  required  to  produce 

the  desired  thermal  efficiency  when  the  coarse  mesh  cores  were  utilized. 

3.  The  coarse  mesh  cores  add  considerably  more  weight  and 
volume  to  the  gas  turbine  arrangement  than  the  fine  mesh  cores  for  a 
given  thermal  efficiency. 

4.  Figure  12  shows  that  required  core  weight  decreases  or 
remains  essentially  constant  as  frontal  area  is  increased  with  7ft A-  ■  «34 
for  all  cores  except  core  4.   It  is  apparent  that  this  should  be  the 
case  after  inspection  of  Figure  11  which  reveals  that  required  volume 
remains  essentially  constant  as  frontal  area  is  increased  for  all  cores 
except  core  4. 

5.  Total  relative  pressure  drop  vs.  regenerator  effective- 
ness appears  to  be  a  linear  relationship  for  fixed  thermal  efficiency 
independent  of  internal  core  geometry. 

6.  Fine  mesh  cores  are  much  more  sensitive  to  changes  in 
external  dimensions  than  coarse  mesh  cores.   Shorter  length  increments 
than  .1  ft  are  necessary  for  good  accuracy  when  investigating  fine  mesh 
cores.  This  is  apparent  from  Figure  13  where  the  points  determined  for 
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cores  3  and  6  are  the  most  seriously  in  error.   Instructions  for  changing 
the  length  increments  in  the  computer  program  are  given  in  Appendix  A. 
7.   Higher  maximum  thermal  efficiencies  are  obtainable 
with  the  coarse  mesh  cores,  however,  the  corresponding  lengths  and  weights 
are  much  greater  than  with  the  fine  mesh  cores. 

8.  Cores  3  and  6,  having  approximately  the  same  plate  spacing, 
produced  essentially  the  same  results  even  though  core  6  had  perforated 
fins  and  3  did  not. 

9.  Core  5  produced  results  similar  to  3  and  6  because  of  its 
extremely  high  friction  factor  even  though  the  plate  spacing  is  nearer 
to  that  of  cores  1,  2  and  4. 

10.  All  cores  produced  a  substantial  increase  in  thermal  effi- 
ciency over  that  obtainable  with  the  open  cycle. 
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FIGURE  11.   SECTIONAL  DRAWING  ILLUSTRATING  CORE  GEOMETRIES 
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Table  I   Characteristics  and  Program  Results  -  Core  1 

Designation,  Plain  plate-fin  surface  10.27  T 

Plate  spacing,  b  s  0.544  in. 

Flow  passage  hydraulic  diameter,  4r^  »  .01259  ft. 

Compactness,  fi    -  289.93  ft2/ft3 

Fin  metal  thickness  ■  .010  in. 

Fin  area/total  area  ■  0.863 

Fin  spacing  at  plates,  c  ■  0.1946  in. 

Results  at  thermal  efficiency  s  .34 


Frontal 

Length 

Regenerator 

Area 

(ft) 

Effectiveness 

3.0 

desired 

thermal  effici 

3.5 

6.2 

.793 

4.0 

4.4 

.754 

4.5 

3.7 

.738 

5.0 

3.3 

.731 

5.5 

2.9 

.720 

6.0 

2.7 

.719 

6.5 

2.5 

.716 

7.0 

2.3 

.710 

vJ^ight 
(lb) 


1122 
910 
816 
853 
825 
838 
840 
832 


AP 


.116 
.079 
.062 
.053 
.046 
.042 
.038 
.036 


total 


Results  at  maximum  thermal  efficiency 


Frontal 

Thermal 

Optimum 

Regenerator 

Weight 

L±r 

Area 

Efficiency 

Length 

(ft) 

Effectiveness 

(lb) 

V  p ) 

3.0 

.320 

5.9 

.765 

985 

.143 

3.5 

.343 

7.3 

.816 

1321 

.139 

4.0 

.359 

8.5 

.850 

1758 

.132 

4.5 

.371 

9.7 

.875 

2257 

.129 

5.0 

.382 

9.7 

.885 

2508 

.114 

5.5 

.391 

10.2 

.897 

2901 

.106 

6.0 

.399 

11.0 

.910 

3413 

.103 

6.5 

.406 

11.7 

.919 

3932 

.100 

7.0 

.411 

12.3 

.927 

4452 

.097 

total 
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Table  II   Characteristics  and  Program  Results  -  Core  2 


Designation,  Plain  plate-fin  surface  12.00  T 

Plate  spacing,  b  =  0.250  in. 

Flow  passage  hydraulic  diameter,  4r, 

Compactness,  /£   -  392.7  ft2/ft3 

Fin  metal  thickness  ■  .006  in. 

Fin  area/total  area  ■  .773 

Fin  spacing  at  plates,  c  ■  .1667  in. 


.009412  ft. 


Results  at  thermal  efficiency  ■  .34 


Frontal 

Length 

Area 

(ft) 

3.0 

desir 

3.5 

3.1 

4.0 

2.6 

4.5 

2.2 

5.0 

2.0 

5.5 

1.8 

6.0 

1.7 

6.5 

1.6 

7.0 

1.5 

Regenerator 
Effectiveness 


Weight 
(lb) 


desired  thermal  efficiency  not  obtainable 

.769  607 

.749  582 

.730  554 

.726  560 

.718  554 

.719  571 

.719  582 

.718  588 


AR 


095 
071 
057 
,049 
044 
040 
037 
035 


total 


Results  at  maximum  thermal  efficiency 


rontal 

Thermal 

Optii 

Area 

Efficiency 

Lengt 

3.0 

.332 

4.0 

3.5 

.349 

4.9 

4.0 

.365 

5.8 

4.5 

.377 

6.6 

5.0 

.388 

7.2 

5.5 

.397 

7.8 

6.0 

.404 

8.3 

6.5 

.410 

8.3 

7.0 

.415 

8.4 

i      Regenerator 
Length  (ft)   Effectiveness 


.803 
.838 
.866 
.887 
.902 
.914 
.924 
.928 
.933 


Weight 
(lb) 

H^ 

672 

.149 

960 

.142 

1299 

.135 

1663 

.128 

2016 

.121 

2402 

.116 

2788 

.111 

2984 

.101 

3292 

.096 

'total 
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Table  III   Characteristics  and  Program  Results  -  Core  3 


Designation,  Plain  plate-fin  surface 

Plate  spacing,  b  -  0.100  in, 

Flow  passage  hydraulic  diair.et.er,  4r^ 

Compactness,  &     -  1332.45  ft2/ft3 

Fin  metal  thickness  -  .002  in. 

Fin  area/total  area  ■  .837 

Fin  spacing  at  plates,  c  ■  .0431 


46.45  T 

-  0.002643  ft 


Results  at  thermal  efficiency  ■  .34 


Frontal 

Length 

Area 

(ft) 

3.0 

desi 

3.5 

desi 

4.0 

0.5 

4,5 

0.3 

5.0 

0.3 

5.5 

0.2 

6.0 

0.2 

6.5 

0.2 

7.0 

0.2 

Regenerator 
Effectiveness 


Weight 
(lb) 


desired  thermal  efficiency  not  obtainable 
desired  thermal  efficiency  not  obtainable 
.819  206 

.760  139 

.779  154 

.733  113 

.750  123 

.765  134 

.778  144 


A_P 

P 


.139 
.084 
.077 
.055 
.051 
.048 
.045 


total 


Results  at  maximum  thermal  efficiency 


Frontal 

Thermal 

Opti 

Area 

Efficiency 

Lengt 

3.0 

.303 

0.5 

3.5 

.323 

0„6 

4.0 

.342 

0.6 

4.5 

.355 

0.7 

5.0 

.368 

0.7 

5.5 

.377 

0.8 

6.0 

.389 

0.7 

6.5 

.397 

0.8 

7.0 

.403 

0.9 

i      Regenerator    Weight 
Length  (ft)   Effectiveness    (lb) 


.769  154 

.820  216 

.841  247 

.872  324 

.884  360 

.905  452 

.904  432 

.920  535 

.933  648 


P  / 


.186 
.185 
.162 
.164 
.148 
.150 
.121 
.124 
.126 
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Table  IV   Characteristics  and  Program  Results  -  Core  4 


Designation,  Plain  plate-fin  surface  6.2 

Plate  spacing,  b  =  .405  in. 

Flow  passage  hydraulic  diameter. 

Compactness,  &     -  204.0  ft2/ft3 

Fin  metal  thickness  =  0.010  in. 

Fin  area/total  area  ■  .728 

Fin  spacing  at  plates,  c  ■  .322  in. 


4rh  -  0.01.82  ft 


Results  at  thermal  efficiency  =  .34 


Frontal 

Length 

Regenerator 

Area 

(ft) 

Effectiveness 

3.0 

6.2 

.756 

3.5 

5.8 

.736 

4.0 

5.7 

.726 

4.5 

5.6 

.721 

5.0 

5.4 

.716 

5.5 

5.2 

.714 

6.0 

4.9 

.712 

6.5 

4.6 

.710 

7.0 

4.3 

.708 

Weight 
(lb) 


887 
968 
1087 
1202 
1287 
1364 
1402 
1426 
1435 


'aP 
P 


.082 
.062 
.052 
.046 
.042 
.038 
.035 
.033 
.031 


total 


Frontal 
Area 


Results  at  maximum  thermal  efficiency 


Thermal      Optimum      Regenerator    Weight 
Efficiency   Length  (ft)  Effectiveness    (lb) 


V  P 


total 


3.0 

.355 

10.9 

.848 

1559 

.137 

3.5 

.371 

13.3 

.869 

2219 

.123 

4.0 

.382 

16.2 

.886 

3090 

.116 

4.5 

.391 

18.5 

.897 

3969 

.106 

5.0 

,399 

20.7 

.907 

4935 

.099 

5.5 

.405 

21.4 

.912 

5612 

.091 

6.0 

.411 

22.4 

.920 

6551 

.086 

6.5 

.415 

22.7 

.923 

7035 

.079 

7.0 

.420 

23.5 

.929 

7843 

.075 
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Table  V   Characteristics  and  Program  Results  -  Core  5 


Designation,  Wavy-fin  plate-fin  surface  11.5-3/8W 

Plate  spacing,  b  =  .375  in. 

Flow  passage  hydraulic  diameter,  4r^  =  0.00993  ft. 

Compactness,  /§    =  345  tt2/ft3 

Fin  metal  thickness  -  .010  in 

Fin  area/total  area  =  .822 

Fin  spacing  at  plates,  c  =  .087  in. 


Results  at  thermal  efficiency  =  .34 


Frontal     Length      Regenerator 
Area        (ft)      Effectiveness 


3.0 
3.5 
4.0 
4.5 
5.0 
5.5 
6.0 
6.5 
7.0 


Weight 
(lb) 


k 


P 


P 


total 


desired  thermal  efficiency  not  obtainable 

desired  thermal  efficiency  not  obtainable 

desired  thermal  efficiency  not  obtainable 

0.9         .772  335  .092 

0.8         .758  331  .073 

0.7         ,739  319  .060 

0.7  .745  347  .054 

0.6         .721  323  .046 

0.6         .726  347  .043 


Results  at  maximum  thermal  efficiency 


Frontal 
Area 


Thermal 
Efficiency 


Optimum 
Length  (ft) 


Regenerator    Weight 
Effectiveness   (lb) 


total 


3.0 

.301 

0.9 

3.5 

.321 

1.1 

4.0 

.338 

1.3 

4.5 

.352 

1.6 

5.0 

.364 

1.7 

5.5 

.374 

2,0 

6.0 

.383 

2.2 

6.5 

.391 

2.3 

7.0 

.398 

2.5 

742 

223 

.168 

788 

319 

.161 

821 

430 

.152 

854 

596 

,153 

866 

703 

.137 

888 

910 

.131 

900 

1092 

.131 

907 

1237 

.120 

916 

1448 

.116 

*Weights  are  low  since  they  were  based  on  core  length  not  flow 
passage  length 
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Table  VI   Characteristics  and  Program  Results  -  Core  6 


Designation,  Perforated-fin,  plate-fin  surface,  Solar  #2 

Plate  spacing,  b  =  0.100  in. 

Flow  passage  hydraulic  diameter,  4r,  ■  .0046272  ft. 

Compactness,  /&   -  673.143  ft2/ft3 

Fin  metal  thickness  =  .0047  in. 

Fin  area/total  area  ■  .645 

Fin  spacing  at  plates,  c  ■  .075  in. 


Results  at  thermal  efficiency  =  .34 


Frontal     Length      Regenerator 
Area        (ft)      Effectiveness 


Weight 
(lb) 


AP 


total 


3.0 
3.5 
4.0 
4.5 
5.0 
5.5 
6.0 
6.5 
7.0 


desired  thermal  efficiency  not  obtainable 

desired  thermal  efficiency  not  obtainable 

desired  thermal  efficiency  not  obtainable 

desired  thermal  efficiency  not  obtainable 

0.7         .768  523        .089 

0.6         .763  493        .074 

0.5          .749  449        .061 

0.4         .725  389       .051 

0.4         .741  419       .048 


Results  at  maximum  thermal  efficiency 


rontal 

Thermal 

Optii 

Area 

Efficiency 

Lengt 

3.0 

.287 

0.8 

3.5 

.305 

1.0 

4.0 

.322 

1.1 

4.5 

.338 

1.2 

5.0 

.352 

1.2 

5.5 

.365 

1.4 

6.0 

.375 

1.5 

6.5 

.386 

1.5 

7.0 

.393 

1.6 

Regenerator    Weight 
Effectiveness    (lb) 


722  359 

766  523 

793  658 

822  808 

840  897 

871  1151 

889  1346 

901  14  53 

913  1675 


aP 
.  P 


total 


186 
179 
164 
153 
136 
140 
136 
125 
123 
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APPENDIX  A 
Computer  Program  Usage 

The  computer  program  is  written  in  FORTRAN  IV  language  described  in 
(9).   The  main  program  and  all  subroutines  were  written  using  double 
precision  in  order  to  obtain  the  best  possible  accuracy.   This  is  parti- 
cularly essential  in  subroutines  where  least  squares  curve  fitting  is 
utilized  to  determine  a  functional  relationship  from  a  set  of  data  points, 

Required  computer  storage  capacity  is  small  as  no  subscripted  vari- 
ables or  arrays  were  introduced  while  writing  the  program.  Essentially, 
the  only  storage  requirement  is  for  program  compilation. 

Approximately  three  minutes  of  computer  time  on  an  IBM  System/360 
is  required  for  each  output  type  described  in  section  six  with  the  excep- 
tion of  output  types  one  and  six  which  require  only  about  30  seconds. 

a.  Subroutine  XMAIN 

The  main  program  is  designed  for  sizing  counterflow  heat 
exchangers  only.   Subroutine  XMAIN  is  general  however  and  could  be  uti- 
lized for  numerous  other  flow  arrangements  provided  that  function  GENFX 
(TU,  RATIO),  called  by  subroutine  XMAIN,  is  changed  to  a  regenerator  ef- 
fectiveness function  corresponding  to  the  desired  flow  arrangement. 

b.  Subroutines  Dependent  on  Internal  Core  Geometry 

The  relationships  of  friction  factor  and  Colburn  j -factor 
vs.  Reynolds  number  must  correspond  to  the  core  geometry  being  investi- 
gated.  These  relationships  are  fed  into  the  program  as  function  sub- 
programs FRICC(RE),  FRICH(RE),  STPRC(RE)  AND  STPRH(RE) .   An  example  is 
shown  on  page  111.   Only  the  cold  side  functions  are  shown  as  the  same 
core  was  used  for  both  hot  and  cold  sides  in  all  cases. 

Subroutine  KEKC  (SIGMA,  RE,  XKE,  XKC) ,  K  and  K  as  a  function  of 
(J~     and  Re,  depends  on  the  shape  of  the  core  flow  passages.   Subroutines 
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were  written  for  square  and  triangular  passages  and  are  shown  on   pages 
105  and  106. 

c.   Data  Cards 

All  general  inputs  for  output  types  one  through  five 
are  read  into  the  computer  program  from  the  first  four  data  cards.   Ten 
columns  are  allotted  for  each  input.   Each  of  the  general  inputs  must 
be  of  floating  point  form.   The  inputs  with  corresponding  units  are 
listed  below  in  the  proper  order: 

compressor  inlet  temperature  °R 

-2 
compressor  inlet  pressure  Ibf-ft 

compressor  pressure  ratio 

compressor  efficiency 

turbine  inlet  temperature  °R 

turbine  efficiency 

assumed  regenerator  effectiveness 

fuel  inlet  temperature  °R 

combustor  efficiency 

lower  heating  value  of  fuel  BTU-lbm 

-2 

hot  side  regenerator  exhaust  pressure   Ibf-ft 

horsepower  output  hp 

assumed  relative  pressure  drop 
through  regenerator  cold  side 

assumed  relative  pressure  drop 
through  regenerator  hot  side 

gas  constant  for  air  ft-lbf (lbm-°R) 

regenerator  plate  thickness  in 

cold  side  transfer  area/volume  1 

between  cold  side  plates  ft 
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hot  side  transfer  area/volume  ft-lbf (lbm-°R) 

between  hot  side  plates 

cold  side  plate  spacing  in 

hot  side  plate  spacing  in 

cold  side  hydraulic  radius  ft 

hot  side  hydraulic  radius  ft 

mass  flow  rate  lbm-hr 

relative  pressure  drop  through 
combustor 

cold  side  fin  thickness  ft 

hot  side  fin  thickness  ft 

2       -1 
fin  thermal  conductance  BTU-(hr-ft  -°F/ft) 

cold  side  fin  or  extended  area/ 
total  area 

hot  side  fin  or  extended  area/ 
total  area 

-3 
density  of  regenerator  material        lbm-ft 

If  output  type  six  is  desired,  only  the  first  12  general  inputs  are 

necessary.   The  third  and  fourth  data  cards  must  be  placed  in  the  deck$ 

however,  they  may  be  left  completely  blank.   The  data  cards  utilized  in 

generating  the  results  for  core  2  are  shown  on  page  95  . 

The  fifth  data  card  is  used  to  designate  the  type  of  output  desired. 
The  fixed  point  integer  "1"  in  column  one  designates  output  type  one; 
the  fixed  point  integer  "2"  in  column  one  designates  output  type  two 
etc. 

The  sixth  and  last  data  card  depends  on  the  type  of  output  designat- 
ed.  In  all  cases  numbers  must  be  in  floating  point  format.   Ten  columns 

are  allowed  for  each  input  beginning  in  column  one.   Units  are  either  ft. 

2 
or  ft.   corresponding  to  length  or  area.   The  inputs  for  each  type  out- 
put are  listed  below  in  the  proper  order: 
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output  type  one:         length 

frontal  area 
output  type  two:         shortest  length 

longest  length 
length  increment  size 
thermal  efficiency 
output  type  three:       smallest  area 

largest  area 
area  increment  size 
thermal  efficiency 
output  type  four:        thermal  efficiency 
output  type  five:        smallest  area 

largest  area 
area  increment  size 
output  type  six:         no  data  card 
d.   Incidental  Instruction  in  Program  Usage 

There  are  several  potential  program  trouble  areas  that 
the  user  should  be  aware  of  when  investigating  various  core  geometries. 
First  of  all,  it  is  important  to  keep  in  mind  the  range  of  Reynolds 
numbers  over  which  the  functional  relationships  of  friction  factor  and 
Colburn-j  factor  have  been  obtained.   Frontal  areas  must  be  chosen  so  that 
the  hot  and  cold  Reynolds  numbers  fall  in  this  range.   This  is  of  course 
largely  dependent  on  internal  regenerator  core  geometry  and  mass  flow 
rate.   Secondly,  experience  has  shown  that  the  iteration  process  diverges 
when  the  total  relative  pressure  drop  through  both  sides  of  the  core  and 
combustor  exceeds  approximately  307o.   This  pressure  drop  is  mainly  a 
function  of  regenerator  length  and  core  plate  spacing.   In  general,  shorter 


89 


lengths  are  necessary  as  plate  spacing  is  decreased. 

Main  program  statements  061,  125,  126  and  141  designate  the  initial 
frontal  area  and  length  used  in  incrementing  out  the  desired  values  for 
output  types  two  and  four.   These  statements,  as  shown  in  the  program 
printout,  are  set  for  a  600  hp  gas  turbine  with  a  regenerator  having  a 
core  plate  spacing  of  .2  inches  or  greater.   For  cores  with  plate  spacing 
smaller  than  .2  inches  or  exceptionally  high  friction  factor  it  may  be 
necessary  to  shorten  the  length  in  statement  125  and  to  increase  the 
frontal  area  in  statements  061,  126  and  141  to  get  the  incrementing  pro- 
cess started  properly  It  may  also  be  necessary  to  increase  the  frontal 
area  in  the  above  statements  for  larger  gas  turbine  arrangements. 

Main  program  statements  62,  92,  125,  140,  153  and  177  designate  the 

length  and  area  increments  to  be  used  in  the  main  program  output  types. 

2 
They  are  presently  set  for  .1  ft  and  .1  ft  but  may  be  changed  as  neces- 
sary if  better  accuracy  is  desired. 

e.   Listing  of  Program  Variables 

Al       dummy  variable  used  in  first  read  statement 

A2       dummy  variable  used  in  first  read  statement 

AFINC     cold  side  fin  area/total  area 

AFINH     hot  side  fin  area/total  area 

AFLOWC    cold  side  free-flow  area 

AFLOWH    hot  side  free -flow  area 

AFRINC    incremental  frontal  area 

AFRON     frontal  area 

AFR0N1    smallest  frontal  area  considered 

AFR0N2    largest  frontal  area  considered 

AHTC      cold  side  heat  transfer  area 

AHTH      hot  side  heat  transfer  area 
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ALFAC 

ALFAH 

BC 

BETAC 

BET  AH 

BH 

CAPC 

CAPH 

CEFF 

COBJC 

COBJH 

COMEFF 

CPC 

CPH 

CRATIO 

DEN 

FC 

FCOND 

FH 

FINEFC 

FINEFH 

FINTHC 

FINTHH 

GC 

GENEFF 

GH 

HI 


total  transfer  area  of  cold  side  of 
regenerator/total  regenerator  volume 

total  transfer  area  of  hot  side  of 
regenerator/total  regenerator  volume 

cold  side  plate  spacing 

cold  side  transfer  area/cold  side  volume 

hot  side  transfer  area/hot  side  volume 

hot  side  plate  spacing 

cold  side  capacity  rate 

hot  side  capacity  rate 

compressor  efficiency 

cold  side  j  factor 

hot  side  j  factor 

combustor  efficiency 

specific  heat  of  cold  side  air 

specific  heat  of  hot  side  air 

compression  ratio 

density  of  core  material 

cold  side  friction  factor 

fin  thermal  conductance 

hot  side  friction  factor 

cold  side  fin  efficiency 

hot  side  fin  efficiency 

cold  side  fin  thickness 

hot  side  fin  thickness 

cold  side  mass  flow  rate  per  unit  area 

regenerator  effectiveness 

hot  side  mass  flow  rate  per  unit  area 

enthalpy  of  air  at  compressor  inlet 
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H2 

H3 

H4 

H5 

H2S 

H5S 

HEATC 

HEATH 

HETVAL 

HF 

HP 

PI 

P2 

P3 

P4 

P5 

P6 

P23 

P56 

P23R 

P34R 
P56R 

PR1 

PR2S 

PR5S 

PR4 

PRC 


compressor  discharge  enthalpy 

combustor  inlet  enthalpy 

turbine  inlet  enthalpy 

turbine  exhaust  enthalpy 

compressor  discharge  ideal  enthalpy 

turbine  exhaust  ideal  enthalpy 

cold  side  heat  transfer  coefficient 

hot  side  heat  transfer  coefficient 

lower  heating  value  of  fuel 

fuel  enthalpy 

horsepower  output 

compressor  inlet  pressure 

compressor  discharge  pressure 

combustor  inlet  pressure 

turbine  inlet  pressure 

turbine  exhaust  pressure 

regenerator  hot  side  exhaust  pressure 

pressure  drop  through  regenerator  cold  side 

pressure  drop  through  regenerator  hot  side 

relative  pressure  drop  through  regenerator 
cold  side 

relative  pressure  drop  through  combustor 

relative  pressure  drop  through  regenerator 
hot  side 

relative  pressure  at  compressor  inlet 

relative  pressure  at  point  2s. 

relative  pressure  at  point  5s. 

relative  pressure  at  turbine  inlet 

cold  side  Prandtl  number 
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PRH 

R 

RADC 

RADH 

RATIO 

REC 

REH 

RELP 

RFA 

SFC 

SIGMAC 

SIGMAH 

STC 

STH 

SUREFC 

SUREFH 

SVOL 

Tl 

T2 

T3 

T4 

T5 

T6 

T2S 

T5S 

TCM 

TEFF 


hot  side  Prandtl  number 

gas  constant  for  air 

cold  side  hydraulic  radius 

hot  side  hydraulic  radius 

capacity  rate  ratio 

cold  side  Reynolds  number 

hot  side  Reynolds  number 

total  relative  pressure  drop 

fuel  to  air  ratio 

specific  fuel  consumption 

cold  side  free  flow  area/cold  side 
frontal  area 

hot  side  free  flow  area/hot  side 
frontal  area 

cold  side  Stanton  number 

hot  side  Stanton  number 

cold  side  surface  efficiency 

hot  side  surface  efficiency 

solid  core  volume 

compressor  inlet  temperature 

compressor  discharge  temperature 

combustor  inlet  temperature 

turbine  inlet  temperature 

turbine  exhaust  temperature 

regenerator  hot  side  exhaust  temperature 

temperature  at  point  T_ 

temperature  at  point  T_ 

mean  cold  side  temperature 

turbine  efficiency 
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TF 

THEFF 

THM 

TU 

U 

V2 

V3 

V5 

V6 

VCM 

VHM 

VIC 

VIH 

VOL 

W 

WT 

XKCC 

XKCH 

XKEC 

XKEH 

XL 

XL1 

XL2 

XLINC 

XMDOT 


fuel  temperature 
thermal  efficiency 
mean  hot  side  temperature 
number  of  heat  transfer  units 
overall  heat  transfer  coefficient 
cold  side  inlet  specific  volume 
cold  side  outlet  specific  volume 
hot  side  inlet  specific  volume 
hot  side  outlet  specific  volume 
mean  cold  side  specific  volume 
mean  hot  side  specific  volume 
cold  side  viscosity 
hot  side  viscosity 
regenerator  volume 
regenerator  plate  thickness 
regenerator  weight 
cold  side  contraction  coefficient 
hot  side  contraction  coefficient 
cold  side  expansion  coefficient 
hot  side  expansion  coefficient 
regenerator  length 
shortest  length  considered 
longest  length  considered 
length  increment 
mass  flow  rate 
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APPENDIX  B 
Curve  Fitting 

It  is  frequently  necessary  to  represent  by  a  functional  relation- 
ship data  that  is  given  in  the  form  of  a  set  of  x-y  points.   One  may 
want  to  know  the  value  of  "yM  corresponding  to  some  "x"  that  lies  between 
two  tabulated  "x"  values  or  may  want  to  find  the  "y"  value  corresponding 
to  an  "x"  value  outside  the  range  of  the  tabulated  "x"  values.  The 
problem  then  is  to  find  a  curve  that  approximates  the  given  data  with 
sufficient  accuracy. 

Before  the  computer  program  could  be  fully  developed,  various  func- 
tional relationships  of  the  type  mentioned  above  had  to  be  obtained. 
Many  of  the  necessary  relationships  were  available  in  literature 5  how- 
ever, seven  of  them  were  not  and  had  to  be  obtained  by  some  appropriate 
curve  fitting  method.  The  "least  squares"  method  was  utilized  in  fitting 
six  of  the  curves  and  straight  line  approximations  were  used  on  the  re- 
maining one.   It  was  decided  to  use  the  least  squares  method  because  of 
its  simplicity  and  the  smoothness  of  fit  that  is  realizable  in  most 
cases. 

a.  Least  Squares  Theory. 

The  deviation  at  a  data  point  is  defined  as  the  difference 
between  the  tabulated  y-value  and  the  y-value  computed  from  a  functional 

relationship.  The  least  squares  method  asks  for  a  minimum  of 

N 

"V  /      .2- 

(B-l) 


Yf'j'c-'fsj 


where  N  equals  the  number  of  data  points,  y'  denotes  a  value  computed 
from  the  functional  relationship  and  y  denotes  the  actual  tabulated  value 
This  leads  to  equations  that  in  many  cases  of  practical  interest  are 
linear  and,  in  principle  at  least  easy  to  solve,  If  it  is  assumed  that 
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the  functional  relationship  desired  is  of  the  form 


(B-2) 


a  ,  a,,  a.  a.  must  be  determined  to  minimize 

o   1   I  k    ki 

2- 


S=£(*(XtfW0 


(B-3) 


(.*■! 


Taking  the  derivative  of  S  with  respect  to  a  and  setting  it  equal  to  zero 


0* 


_    £1-    } 
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r=»  £~ 


a/ 
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"  I  iJL  iA 


1 £.1     / 
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I * — '         /  4—  fB' 


4) 


This  gives  one  equation.   By  taking  the  derivatives  of  S  with  respect  to 
the  other  coefficients,  a.,  through  a,  »  another  "k"  equations  are 


V-i 


(-/ 


<=/ 


(B-5) 


/V    ,     ,  A/     .  ,    >J     'L& ,  M 
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The  k+1  equations  can  then  be  solved  for  the  k+1  coefficients. 

The  value  selected  for  k  can  be  any  value  less  than  N.   For  best 
results,  however,  k  should  be  made  only  large  enough  to  give  desired 
accuracy.   As  k  is  increased  the  computer  round-off  errors  inherent  in 
solving  the  set  of  simultaneous  equations  becomes  more  significant  and 
accuracy  is  lost  not  gained.   The  generated  curve  may  also  tend  to  fluc- 
tuate between  data  points  if  k  is  too  large,  making  a  smooth  fit  impos- 
sible. 

b.   Gauss  Elimination. 

The  method  of  Gauss  elimination  was  used  for  solving  the 
simultaneous  equations  mentioned  above.   A  system  of  three  equations  in 
three  unknowns  will  be  used  to  illustrate  the  procedure  because  the 
nomenclature  becomes  somewhat  complicated  for  the  general  case  of  "n" 
equations  and  "n"  unknowns.   Let  x. ,  x  and  x.  be  the  unknowns  and  the 
equations  are 

(B-6) 

The  equations  must  be  arranged  so  that  a. 1  is  not  zero.   Next  define  a 
multiplier, 


AAA  -   ^l-' 


All  <B"7> 


The  first  equation  is  then  multiplied  by  m_,  and  subtracted  from  the 
second  equation.   The  result  is 
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But 

so  x.  has  been  eliminated  from  the  second  equation. 

Define: 

/ 
tfii.  =  #n.-**,.ai2_  (B-9) 

/ 

#*3  =  fly  -/»*»!  a»3  (B-10) 

A'  =  A  -  /**,/,  (B-ll) 

then  equation  (B-8)  becomes 

If  a  multiplier  is  now  defined  for  the  third  equation, 

^  -  4ft.  (,.13) 

and  the   first  equation  is   subtracted   from  the   third,    the  result  is 

#3i'*z.  **jj'Xj  *  ^  (B-14) 


where 


C\ir    -  ^3u-"  />*\j  #,x  (B-15) 


^32'  =  tfjj  —  #n3a  13  (B-16) 

^   «  /£3  -  /M3  /,  (B-17) 
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The  resulting  three  equations  in  three  unknowns  are 

flu*,  -*-a,i.Xi.+tf  ,3X3  =  4-t 

*n'xt.+a*l*j-ti  (B"18) 

AixXx  +a33/xJ-  li 

If  X.  is  now  eliminated  from  the  last  equation  using  the  same  procedure 
the  following  is  obtained: 

Aw  X|  +  aivXj.+-^ijX3  -  Jtx 

AvL-Ki  +  txi'li  -  l!  (B-19) 

&*{  X3  ^  1/ 

It  is  now  a  straight  forward  process  to  solve  for  x_,  to  substitute  that 
result  into  the  second  equation  to  find  x_,  and  finally  to  substitute  in- 
to the  first  equation  to  find  x-  . 

#33 
X-   (  ""*  ~^lZ   Xj'  (B-20) 

&XX, 

3      An 


c.  Curve  Fitting  Computer  Program 

The  computer  program  shown  on  the  following  pages  was 
utilized  in  fitting  the  necessary  curves.  The  main  program  manipulates 
the  data  point  information  into  the  equation  form  shown  in  part  "a"  of 
this  appendix.   Subroutine  GAUSS  2  is  then  called  and  the  equations  are 
solved  simultaneously  using  the  method  described  above  (10).   The  coeffi' 

cients,  a, a,  ,  determined  for  a  polynomial  of  the  form  y  «  a.  + 

2  k-1 

a„x  +  a~x  a.x    are  printed  out  along  with  a  column  of  x-y 
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points  generated  from  the  polynomial.  These  points  can  then  he  compared 
with  the  original  data  points  to  determine  whether  or  not  the  fit  is 
sufficiently  accurate. 

A  maximum  of  fifty  data  points  can  be  read  into  the  computer.  For 
best  results,  they  should  be  as  evenly  spaced  along  the  "x"  axis  as  pos- 
sible. A  special  problem  arises  when  using  the  program  to  fit  Golburn-j 
and  friction  factor  vs.  Reynolds  number  data.  This  data  is  normally 
plotted  on  logarithmic  coordinates  and  is  therefore  not  evenly  spaced  on 
linear  coordinates.  To  remedy  this  situation,  the  natural  logarithms  of 
the  "x"  and  "y"  values  can  be  taken  before  starting  the  curve  fitting 

procedure.   The  resulting  polynomial  is  of  the  form  ln(y)  -  a.  +  aJln  x)  + 

2  k-1 

a~(ln  x)   ai,(-'-n  x)    •   Anti- logarithms  can  then  be  taken  to  deter- 
mine the  desired  "y"  value.   These  additional  steps  are  incorporated  into 
the  computer  program  as  it  is  shown.   If  the  user  desires  to  fit  data 
plotted  on  linear  coordinates,  cards  150,  160,  510,  and  520  must  be  re- 
moved and  card  530  must  be  changed  to  read: 
WRITE  (6,55)  Xl(LQ).  Yl(LQ) 

Coefficients  shown  in  parts  1  and  2  of  Table  VII  are  for  polynomials 
of  the  first  form  shown  above.   The  remaining  coefficients  in  Table  VII 
are  for  polynomials  of  the  form  just  described. 

When  utilizing  the  program,  the  following  data  cards  are  necessary: 
(1)  The  number  of  data  points  to  be  used  in  fitting  the 
curve  is  punched  into  the  first  three  columns  on  the  first  data  card. 
The  order  of  the  polynomial  desired  is  punched  into  the  second  three 
columns  of  the  first  data  card.  Both  of  these  numbers  must  be  in  fixed 
point  integer  form.   The  order  of  the  polynomial  can  be  any  number  up  to 
one  less  than  the  number  of  data  points. 
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(2)  The  x-y  values  of  the  data  points  are  punched  into 
the  next  N  data  cards,  one  data  point  per  card.   The  x-value  is  punched 
into  the  first  ten  columns;  the  y-value  is  punched  into  the  second  ten 
columns.   Both  numbers  must  be  of  floating  point  form. 

(3)  As  mentioned  above,  the  printout  includes  a  column 
of  x-y  points  generated  from  the  computed  polynomial.   The  first  point 
is  computed  using  the  smallest  x-value  and  the  last  point  is  computed 
using  the  largest  x-value  of  the  original  data  points.   The  size  of  the 
desired  increments  between  these  two  data  points  is  punched  into  the  next 
data  card  in  the  first  ten  columns.  Floating  point  form  must  be  used. 

(4)  If  the  user  desires  to  fit  more  than  one  set  of 
points  he  can  place  the  next  set  of  data  cards,  for  those  points,  be- 
hind the  set  just  described. 

(5)  The  next  and  last  data  card  must  contain  the  integer 
2  punched  into  column  three.  This  signals  that  all  data  point  sets  have 
been  fit  and  stops  the  program. 

d.   Curve  Fitting  Results 

A  sample  of  least  squares  curve  fitting  is  shown  in  Figure 

18.  All  data  points  used  in  curve  fitting  and  the  resulting  polynomial 
coefficients  are  shown  in  Table  VII. 

The  curve  fit  by  straight  line  approximations  is  shown  in  Figure 

19,  with  the  straight  lines  drawn  in.  Two  straight  lines  were  used,  one 
covering  the  range,. 3  <  (J""  <^  .55  the  other  covering  the  range 

.55<  C7~<^.75.  This  approximation  was  used  rather  than  the  least  squares 
theory  because  small  errors  in  C  cause  virtually  no  error  in  the  comput- 
er program  results.  The  curve  is  from  Kays  (5). 
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Table  VII   Data  Points  Used  in  Curve  Fitting  and  Resultin< 

PolynoiTu  I  Coefficients  Used  in  Program  Subroutines 

1.   Dynamic  Viscosity  of  Air  at  Atmospheric  Pressure  from 

y0  a  L716  k   L0~7  poise 

lynorni&l  Coefficients 

-  4132588012  D-01* 

2  ?06898585  D-02 

-.  1318733093  D-05 

.7538864121  D-12 

7^21455415  D-12 

-.5752187450  D-15 

2117993215  D-18 

-.3862860359  D-22 

,2782105119  D-26 

.6763288488  D-32 


Hil: seii  rath 

T°R 

/  Vo 

180 

0.4038 

540 

1,7  742 

900 

i.556 

1260 

1.7  58 

1620 

2.271 

1980 

2.562 

2340 

2,824 

2  700 

3.066 

3060 

3.290 

3420 

3.501 

A3 
A4 

A5 

a6 

a; 

AH 

A9 
A10 


/..      P"andtl  Number  of  Air  from  Hilsenrath  (6) 


T°R 

Prandtl  No 

180 

0,770 

360 

0.739 

540 

0.708 

720 

0.689 

900 

0.680 

1080 

0.680 

12  60 

0.684 

1440 

0.689 

1620 

0.696 

1800 

0.702 

Po lynomial  C oe  £ ficients 


Al 
A2 
A3 

A4 
A5 
A6 
A7 
A8 
A9 
AiO 


7865630968 
2871584338 
3666074673 
1373281188 

15699040?6 
1973361281 
1041204710 
2047252924 
1538501631 


D-00 

D-Qu 
D-06 
D-09 
D-ll 
D-14 
D-17 
D-21 
D-25 


79605304  27  D-29 


Core-  1  (Plain  plate-fin  surface  10.27  I)  Friction  Factor 
from  Kays  and  London  ( 


Reynolds 

Friction 

Number 

Factor 

300 

.066 

400 

.051 

500 

.042 

700 

.032 

900 

,026 

1250 

,0185 

1500 

0165 

2000 

.0135 

3000 

.011 

4000 

.0095 

6000 

.0082 

J v  n om La 1  Coe ff :c  ients 

Al  ■  ,1060442201  D+03 

A2  «  -. 7242808930  D+02 

A3  -  .1760821881  D+02 

A4  =  -.11  71044519  D+01 

A5  -  -.2766351384  D+00 

A6  ---  .6028901137  D-01 

A?  -  -. 1810176966  D-02 

A8  -  -  7  7  28  7807  35  D-03 

A9  *  .1143227913  D-0'1 

AiO  *  -.640920187  3  D-05 

All  ■  . 1272474507  D-06 


*A1  =  -.4132588012  x  1C*1 
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Core  1  (Plain  plate-fin  surface  10.27  T)  Coburn-j 
Factor  from  Kays  and  London  (1) 


Reynolds 

Coburn-j 

Number 

Factor 

300 

.016 

400 

.013 

500 

.0115 

700 

.0037 

900 

.0072 

1250 

.006 

1500 

.00515 

2000 

.0045 

3000 

.00375 

4000 

.00345 

6000 

.0032 

Polynomial  Coefficie n t s 


Al 

a 

-.9478036457 

D+01 

A2 

= 

.5133725600 

D+01 

A3 

= 

-.1236662988 

D-K)l 

A4 

= 

-.2136908795 

D+00 

A5 

= 

.1534935252 

D+00 

A6 

= 

-.  225t;Q9c;^0 

T)-0] 

A7 

= 

-.3614040571 

D-03 

A8 

= 

.3150813943 

D-03 

A9 

= 

-.1222790555 

D-04 

A10 

= 

-.1473493662 

D-05 

All 

■ 

.9202049786 

D-07 

Core  2  (Plain  plate-fin  surface  12.00  T)  Friction 
Factor  from  Kays  and  London  (1> 


Reynolds 

Friction 

Number 

Factor 

200 

.080 

300 

.057 

500 

.037 

700 

.027 

1000 

.019 

1500 

.015 

2000 

.013 

3000 

.0115 

4000 

.0112 

6000 

.0095 

8000 

.0036 

Polynomial  Coefficients 


Al 

B 

.2855654237 

D+03 

A2 

= 

-.4366648607 

D+02 

A3 

= 

-.7695241634 

D+02 

A4 

= 

.3558930941 

D+02 

A5 

= 

-.5846569710 

D+01 

A6 

= 

.2516782204 

D+00 

A7 

= 

.3534066422 

D-01 

A8 

= 

-.3719622996 

D-02 

A9 

= 

-.3919235358 

D-04 

A10 

= 

.1641370472 

D-04 

All 

= 

-.5035864719 

D-06 

Core  2  (Plain  plate-fin  surface  12.00  T)  Coburn-j  Factor 
from  Kays  and  London  ( 1) 


Polynomial  Coefficients 


Reynolds 

Coburn-j 

Number 

Factor 

200 

.016 

300 

.012 

500 

.008 

700 

.0064 

1000 

.0049 

1500 

.0038 

2000 

.00355 

3000 

.0036 

4000 

.0035 

6000 

.00325 

8000 

.003 

Al 

= 

-.1678025363 

D+04 

A2 

- 

.1147089449 

D+0& 

A3 

= 

-.2302525106 

D+03 

A4 

= 

-.9930779398 

D+01 

A5 

= 

.8790472874 

D+01 

A6 

= 

-.7254951992 

D+00 

A7 

= 

-.5905350854 

D-01 

A8 

= 

.6125863884 

D-02 

A9 

= 

.9700663367 

D-03 

A10 

= 

-.1394667289 

D-03 

All 

= 

.4807945572 

D-05 
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Core  3  (Plain  plate-fin  surface  46,45  T)  Friction  Factor 

from  Kays  and  London   1 ■ 


Reynolds 
Number 

300 

400 

500 

700 

900 

1250 

1500 

2000 

3000 

4000 

6000 


rice,  ion 
Factor 

.059 

.045 

.03  7 

.0265 

.020 

.0163 

.0135 

.0128 

.0115 

.0112 

.0110 


Polynomial  Coefficients 


Al 

= 

5218554  745 

D+02 

A2 

s 

-.3337432034 

EH-02 

A3 

- 

.8259113372 

D+01 

A4 

~ 

-.180226^981 

D+01 

A5 

~ 

.4388704605 

DK)0 

A  6 

- 

-4819876576 

D-01 

A? 

= 

-.3056062581 

D-02 

A8 

as 

.7437300503 

0-03 

A  9 

as 

.3659712646 

D-04 

A10 

- 

-.1134  921023 

D-04 

All 

= 

.5000930552 

D-06 

8.   Core  3  (Plain  place-fin  surface  46,45  T)  Coburn-j  Factor 

from  Kays  and  London  (i ? 


Polynomial  Coefficients 


Key no j  da 

L'oburn-  i 

Number 

Factor 

300 

.0135 

400 

,0105 

500 

.0036 

700 

,0065 

900 

.0053 

12  50 

.0041 

1500 

.0036 

2000 

.00295 

3000 

.0026 

4000 

.0025 

6000 

.00225 

Al 

=3 

-.9179528233 

D+01 

A  2 

=2 

.2549082428 

D+01 

A3 

«* 

-.1138319133 

D+01 

A4 

= 

.1312317457 

D+01 

A5 

as 

-.52007  20093 

D+00 

A6 

a 

.b707637499 

D-01 

A7 

- 

.2230241152 

D-02 

A8 

= 

-.1010311829 

D-02 

A9 

- 

.1206538219 

D-04 

ilO 

- 

.8449625214 

D-05 

ill 

JS 

-.4529482542 

D-06 

9.   Core  4  (Plain  plate-fi-  surface  6.2)  Friction  Factor 
from  Kays  and  London  I 


Pol  y  n  om  la  I  Coe  f  f  i  c  1  e  n_t  s 


Reynolds 

Friction 

Number 

Factor 

.077 

:;0>'^ 

,022 

1250 

.015 

2000 

.010 

2300 

,0098 

2700 

.0093 

3000 

.0092 

4000 

..  008  )  5 

6000 

,0081 

8000 

.00775 

10000 

.0073 

Al 

a 

.4512980279 

D+05 

- 

-.3083071127 

D+05 

A3 

= 

.7727049904 

D-K)4 

A4 

= 

-. 7410141046 

D+03 

A5 

a 

-.6317170678 

D+01 

A6 

= 

.388064 '509 

D+01 

■' 

= 

.6306627426 

D-01 

A8 

=s 

.321948  7556 

D-01 

A9 

= 

-.1290500496 

D-01 

A10 

as 

,1130983538 

D-02 

All 

= 

-.31S9159559 

D-04 
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10.   Core  4   (Plain  plate-fin  surface  6.2)  Cobum-j  Factor 
from  Kays  and  London  (1 


Reynolds 

Coburn- j 

Number 

Factor 

500 

.0072 

800 

.0059 

1250 

.0043 

2000 

.0032 

2300 

.003 

2700 

.0031 

3000 

.0033 

4000 

.0034 

6000 

.0033 

3000 

.00325 

10000 

.00315 

Po] 

.ynomial  Coefficients 

Al 

= 

.7203819966 

D-K)5 

A  2 

= 

-.4396069217 

D+05 

A3 

= 

.1220436223 

D+05 

A4 

= 

-.1162973846 

D+04 

A5 

= 

-.1005649404 

Df-02 

A6 

= 

.6019720585 

D+01 

A7 

= 

.9239757532 

D-01 

A8 

= 

.5306752808 

D-01 

A9 

= 

-.2033142744 

D-01 

l10 

= 

.1762226828 

D-02 

11 

= 

-.4934822700 

D-04 

11.   Core  5  (Wavy-fin  plate-fin  surface  11.5-3/8W)  Friction 
Factor  from  Kays  and  London  (1 


Polynomial  Coefficients 


Reynolds 

Friction 

Number 

Factor 

400 

.14 

500 

.125 

700 

.098 

900 

.083 

1250 

.073 

1500 

.070 

2000 

.0625 

3000 

.0525 

4000 

.0475 

5000 

.043 

7000 

.0375 

Al 

s 

-.1703315790 

D+05 

A  2 

= 

.1290411355 

D+05 

A3 

= 

-.3694551557 

D+04 

A4 

= 

.4359291610 

D+03 

A5 

= 

-.1981584997 

D+01 

A6 

= 

-.3594302044 

D+01 

A7 

= 

.7340606854 

D-01 

A8 

s 

.3091640055 

D-01 

A9 

as 

-.1603543052 

D-02 

A10 

= 

-.7179253410 

D-04 

All 

= 

.5035665432 

D-05 

12.   Core  5  (Wavy-fin  plate-fin  surface  11.5-3/8W)  Cobum-j 
Factor  from  Kays  and  London  (!) 


Polynomial  Coefficients 


Reynolds 

Coburn- i 

Number 

Factor 

400 

.019 

500 

.018 

700 

.0175 

900 

.01625 

1250 

.015 

1500 

.014 

2000 

.013 

3000 

.012 

4000 

.0096 

5000 

.009 

7000 

.0079 

Al 

=3 

.1845757550 

D+05 

A2 

= 

-.1406316423 

D+05 

A3 

SX 

4071246932 

D+04 

A4 

= 

-.4950898847 

D+03 

A5 

= 

.5485777432 

D+01 

A6 

= 

.3828321238 

D+01 

A7 

= 

-.1208492523 

D+00 

A8 

= 

-.2934686653 

D-01 

A9 

= 

.1796872494 

D-02 

A10 

= 

.5114973685 

D-04 

All 

= 

-.4408556059 

D-05 
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13,  Core  6   (Perforated-fin  plate  fin  surface,  Solar 
'  iction  Factor  f.i  1.(8). 


Reynolds 

r iction 

bar 

['actor 

L50 

200 

.  3  1 

300 

73 

500 

.049 

700 

.035 

000 

.  029 

1  2C'r' 

.0225 

1500 

.020 

2000 

.020 

2500 

.020 

3000 

.020 

Pol y n omj a  I  C o e £ f i c i en t s 


.... 

_ 

1900 

-. 

-.3513698265 

D+04 

A3 

s 

10597?  7565 

D+04 

A4 

= 

-.82628 

A5 

SB 

-.2304768472 

D+02 

SB 

.6576122775 

D  r<)  1 

A/ 

a 

-.512507  2877 

A8 

- 

-.7245182215 

D-02 

- 

.3054392764 

D-02 

L0 

= 

-.  1198421133 

D-03 

All 

- 

.  S775033884 

D-07 

14.   Core  6  (Perforated-fin  piate  fin  surface,  Solar  #2) 
Coburn-i  Factor  from  Pucci  et.al.  (8), 


Reynolds 

Coburn-i 

Number 

Factor 

150 

.04 

200 

.03 

300 

.019 

500 

.0115 

700 

.008 

000 

.0065 

1200 

.0054 

1500 

.  0054 

2000 

.0065 

2500 

.0068 

3000 

.0063 

Polynomial  Coefficients 


Al 

— 

-.7987254165 

D+04 

A2 

- 

.0330003623 

D+04 

A3 

- 

-.  1701306104 

D+04 

A4 

as 

.9227812576 

D+02 

A5 

= 

.3735667832 

D+02 

a6 

as 

-.6449659103 

D+01 

A7 

SB 

.1419750320 

D+00 

A8 

= 

.2445179921 

D-01 

A9 

SB 

.2010489258 

D-02 

A10 

SS 

-5001272149 

D-03 

All 

- 

.2063223608 
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